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ABSTRACT 

The  objective  of  this  experiment  was  to  describe  the  high  frequency  acoustic 
backscatter  brom  simulated  arctic  sea  ice  as  a  function  of  the  physical 
properties  of  the  ice.  In  addition,  diffraction  experiments  were  performed  on 
the  edge  of  an  open  lead.  These  edges  are  acoustically  interesting  because 
saline  ice  both  supports  shear  waves  and  is  highly  absorptive.  Experiments 
were  carried  out  on  urea  ice  grown  under  controlled  laboratory  conditions  in 
the  cold  pit  at  the  Cold  Regions  Research  and  Engineering  Laboratory, 
Hanover,  New  Hampshire.  Normal  incidence  sonar  echo  amplitude  data 
were  collected  at  kilo-Hertz  carrier  frequencies  from  the  underside  of  a 
jTOwing  ice  sheet.  Using  a  plane-wave  reverberation  model,  the  velocity  of 
the  acoustic  wave,  attenuation  and  normal  incidence  reflection  coefficient  of 
the  ice  were  measured  at  15, 17, 20, 25  and  190  kHz.  Analysis  showed  that  the 
urea  ice  had  reflection  ooefhdents  of  about  0.4  in  the  tens  of  kHz  range, 
significantly  higher  than  saline  ice.  Attenuation  at  those  h’equendes  was 
about  0.3  dB/cm.  At  higher  frequendes,  the  acoustic  properties  of  the  urea  ice 
begin  to  approach  those  meastired  on  saline  ice.  Fbr  the  diffraction  phase  of 
the  experiment,  a  slot  was  cut  in  in  the  ice  and  the  transducers  were  moved 
laterally  under  the  edge  of  the  lead.  The  diffracted  and  reflected  echo  arrivals 
were  resolved  and  the  amplitude  of  the  diffracted  signal  was  measured.  A 

good  comparison  was  found  between  the  diffracted  arrival  values  and  _ 

Helmholtz-Kirchoff  theory. 


Due  QUALITY  msPEnrigfg^ 


Accesion  For 


NTIS  CRA&I 
one  TAB 

Uriaii::oui,c.od 

Justificatior; 

By . 

Distribution  / 


n 

□ 


Cones 


Oist 


Av;):. 

Sir.'oial 


T 


Table  of  Contents 


page 


Chapter 

1  Introduction 

1 

Chapter 

2  Background 

4 

2.1  Ice  Texture  and  Morphology 

4 

2.2  Acoustic  Theory 

11 

2.2.1  Transmission  of  Power 

11 

2.2.2  Sonar  equation 

17 

2.2.3  Plane  wave  reverberation  model 

20 

2.23.1  Limitations  of  Model 

23 

Chapter 

3  Experimental  Approach 

25 

3.1  Set  up 

25 

3.2  Calibration 

33 

3.2.1  1  /  R  calibration 

33 

3.2.2  Beam  pattern  calibration 

40 

3.2.2.1  Beam  pattern  theory 

43 

3.3  Observations  of  morpholgy 

47 

Chapter 

4  Growth  Experiment 

55 

4.1  Procedure 

55 

4.2  Precautions 

62 

4.3  Observations 

67 

4.4  Analysis 

77 

Chapter 

5  Diffraction  Experiment 

89 

3.1  Di^ction  Experiment  theory 

89 

5.2  Procedure 

94 

5.3  Observations 

109 

Chapter 

6  Error  Analysis 

115 

Chapter 

7  Condusions 

121 

Chapter 

8  Further  Work 

Appendices 

Bibliography 

124 

iv 


Figure 

11 

List  of  Hgures 

c-axis 

page 

5 

Figure 

2.2 

Comparative  morphology 

7 

Figure 

2.3 

Urea  ice  dendritic  interface 

9 

Figure 

Figure 

2.4 

Saline  ice  dendritic  interface 

10 

3.1 

Ice  Growth  rate 

26 

Figure 

3.2 

Temperature  rate 

27 

Figure 

3.3 

Set  up  for  lower  frequencies 

28 

Figure 

3.4 

Schematic  of  voltage  amplifier 

30 

Figure 

3.5 

Set  up  for  190  kHz 

32 

Figure 

3.6  a 

15  kl^  1  /R  pre/post<experiment  calibration 

34 

3.6  b 

17  kHz  1  /R  pre/post-experiment  calibration 

35 

3.6  c 

20  kHz  1  /R  pre/post-experiment  calibration 

36 

3.6  d 

25  kHz  1/R  pre/post-experiment  calibration 

37 

3.6  e 

190  kHz  1/R  dependence 

38 

Figure 

:.7 

Parallel  beam  pat.  calibration 

41 

Figure 

3.8 

Perpendicular  beam  pat.  calibration 

42 

Figure 

3.9 

Urea  ice  Samples  A-C 

51 

Figure 

3.10 

Urea  ice  Samples  D-F 

52 

Figure 

3.11 

Urea  ice  Samples 

53 

Figure 

3.12 

Urea  ice  Samples  I J 

54 

Figure 

4.1  a 

15kHzp.d.f. 

57 

4.1  b 

17  kHz  p.d.f. 

58 

4.1  c 

20kHzp.d.f. 

59 

4.1  d 

25kHzp.d.f. 

60 

4.1  e 

190kHzp.d.f. 

61 

Figure 

4.2 

Growth  experiment  set-ups 

64 

Figure 

4.3  a 

Representative  waveform 

65 

4.3  b 

Multiple  bursts 

65 

4.3  c 

190kHznuU 

66 

Figure 

4.4  a 

190  kHz  reverberation  model  vs.  experiment 

71 

4.4  b 

20  kHz  reverberation  model  vs.  experiment 

72 

4.4  c 

15  kHz  reverberation  model  vs.  e>q)eriment 

73 

4.4  d 

17  kHz  reverberation  model  vs.  e}q>eriment 

74 

4.4  e 

25  kHz  reverberation  model  vs.  e^q^eriment 

75 

Figure 

4.5 

R.C  comparisons  and  attenuation  comparisons 

82 

Hgure 

4.6  a 

15kHz  2nd>  reflections 

85 

4.6  b 

17kHz  2nd*»*  reflections 

86 

4.6  c 

20klb  Tnd-t-  reflections 

87 

4.6  d 

25kHz  2nd4-  reflections 

88 

Figure 

5.1  a 

Diffraction  Experiment  Set  up 

95 

5.1  b 

Lateral  distance  vs  Surface  dtetance 

96 

5.1  c 

Diffraction  Coordinate  System 

97 

Hgure 

5.2 

Travel  time  for  ice  edge 

100 

Hgure 

5.3  a 

Travel  times  for  water  edge 

101 

5.3  b 

Amplitude  for  water  edge 

102 

page 


Figure  5.4  a  Travel  times  for  reflections  103 

5.4  b  Amplitudes  of  reflection  104 

Figure  5.5  Subtraction  of  imder-ice  waveform  107 

Figtire  5.6  Subtraction  of  under-water  wavfonn  108 

Figure  5.7  a  Diffraction  values  lOps  112 

5.7  b  Diffraction  values  A«  50ps  113 

5.7  c  Diffraction  values  A*  200ps  114 


List  of  Tables 


Table 

3 

1/R  calibrations 

36 

Table 

3 

%  error  in  beam  pattern 

47 

Table 

4 

Summary  of  bulk  a>effrcients 

76 

Table 

4 

Digitisation  data  vs  model  R.C. 

78 

Table 

4 

Digitisation  data  vs  model  Attenuation 

84 

Table 

6 

%  error  in  p.di.*s 

114 

Table 

6 

Reflection  coefficient  error 

115 

Appendices 

Appen* 

1 

1  /R  calibration  data 

Aopen< 

2 

Beam  pattern  Calibration  data 

.  •  3pen* 

3 

Plane  wave  reverberation  model 

^  open* 

4 

Diffraction  values 

A.  ^peni 

5 

Diffrraction  theoretical  model 

1 


1.  Introduction: 

Recent  laboratory  experiments  on  acoustic  propagation  through  thin  saline 
ice  (Stanton  et  al.,  1986;  Jezek  et  al.,  1989)  have  demonstrated  that  the  acoustic 
properties  of  saline  ice  are  strongly  correlated  with  the  ice  texture  and 
morphology.  These  experiments  have  shown  that  for  frequencies  above  100  kHz, 
coupling  into  the  ice  is  very  good  (reflection  coefficients  are  typically  less  than 
0.1)  and  the  attenuation  through  the  ice  is  high  (3  dB/cm).  Based  on  statistical 
analysis  of  echo  amplitude  data,  they  conclude  that  these  acoustic  effects  are 
associated  with  the  dendritic  ice-water  interface  of  growing  sea  ice  and  the 
entrapment  of  brine  up  into  the  body  of  the  ice  sheet 

The  sea  ice  cover  of  the  polar  oceans  can  develop  complicated  geometries. 
Scattering  from  the  edges  of  opwn  leads  (chaimels  of  water)  or  from  irregular 
blocks  that  comprise  the  keels  of  sea  ice  ridges  are  important  factors  in  acoustic 
propagation.  Ifigh  frequency  acoustic  backscatter  from  flat  sea  ice  has  been 
studied  since  the  i960's  by  CS.  Qay  and  H.  Medwin  (1977  pp544),  J.R.  Brown  et  al. 
(1966,1967),  G.R.  Garriscm  et  al.  (1978),  and  R.E  Francois  (1988)  and  others.  Only 
recently,  however,  has  theoretical  and  experimental  work  focused  on  backscatter 
from  more  complex  topographical  features  such  as  ice  keels.  A  sophisticated 
model  of  under-ice  acoustic  backscatter  firom  the  Arctic  ice  sheet  has  been 
developed  by  G.C  Bishop  (1989),  and  other  studies  by  Posey  et  al.  (1985)  and  Chin 
Bing  (1985, 1987)  have  reported  under-ice  propagation  and  scattering. 
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The  primary  objective  of  this  experiment  was  to  describe  high  frequency 
acoustic  backscatter  from  simulated  arctic  sea  ice  as  a  function  of  the  physical 
properties  of  the  ice.  Dif^ction  experiments  were  perf(»med  on  the  edge  of  a 
large  hole  cut  in  the  ice  sheet  The  edges  oi  open  leads  are  acoustically  interesting 
because  saline  ice  both  supports  shear  waves  and  is  highly  absorptive.  However, 
to  imderstand  the  diffraction  mechanism  it  was  essential  to  characterize  the  bulk 
acoustic  properties  of  the  ice. 

The  experiment  was  carried  out  tmder  controlled  laboratory  conditions  in 
the  indoor  facilities  at  the  U.S.  Army  Cold  Regions  Research  and  Engineering 
Laboratory,  Hanover,  New  Hampshire.  Urea  ice  is  widely  used  as  a  modelling 
material  for  saline  ice  in  laboratory  tests  because  it  is  non-corrosive  and  has 
similar  mechanical  properties.  In  addition,  the  acoustic  properties  of  ice  e 
believed  to  be  dependent  on  morphology,  and  independent  of  ice  chemistry.  We 
used  urea  ice  as  a  substitute  for  saline  ice  in  our  experiment  since  they  are 
considered  morphologically  very  similar. 

The  experiment  was  divided  into  two  parts: 

a)  The  growth  phase:  Reverberation  experiments  were  performed  on  the 
growing  ice  sheet.  Normal  incidence  sonar  echo  amplitude  data  were  collected  at 
frequencies  of  15kHz,  17  kHz,  20  kHz,  25  kHz  and  190  kHz.  Measurements  were 
taken  every  0.5  cm  of  growth  in  thickness.  The  echo  amplitudes  were  tabulated 
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as  a  function  of  ice  thickness  and  then  fitted  to  a  plane-wave  reverberation 
model.  The  attenuation,  reflection  coefficient  and  phase  velocity  of  the  acoustic 
wave  in  the  ice  were  then  determined. 

b)  The  diffiraction  phase:  Echo  amplitudes  were  measured  as  the  transducers 
were  moved  along  a  horizontal  plane  beneath  the  ice  edge.  The  various  echoes 
from  the  ice  surface,  water  surface,  and  firom  the  ice  edge  were  observed.  The 
received  waveforms  were  then  digitized  and  the  diffracted  arrival  isolated  from 
the  other  signals.  The  amplitude  of  the  diffi'acted  arrival  was  recorded  and  fou^d 
to  be  in  agreement  with  values  predicted  using  Helmholtz-Kirchoff  theory. 
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2.  Background: 

2.1  Ice  Texture  and  Morphology: 

In  recent  years,  urea  ice  has  been  widely  used  for  ntodeliing  sea  ice  in  test 
tanks  in  North  America  and  Europe  (Gow,  1984).  When  urea-doped  water 
freezes,  the  build-up  of  urea  at  the  ice-water  interface  eventually  becomes  too 
large  to  sustain  a  planar  geometry.  The  interface  then  develops  a  ceUular,  or 
doidritic,  morphology.  After  a  stable  dendritic  interface  forms,  further 
freezing  leads  to  the  growth  of  dendrite  tips  which  protrude  down  into  the 
urea  solution.  Simultaneously,  urea  enriched  liquid  becomes  entrapped  in 
the  grooves  between  the  dendrites.  This  mechanical  entrapment  of  urea  is 
mostly  confined  to  the  grooves  between  the  dendrites.  The  dendrites 
themselves  are  elongated  parallel  to  the  crystallographic  basal  plane,  that  is, 
normal  to  the  c-axis.  (See  Figure  2.1)  The  c-axis  is  the  principal  hexagonal  axis 
of  the  ice  crystal;  it  is  widely  used  as  the  reference  axis  for  describing  crystal 
structure  and  related  physical  and  optical  properties  of  the  ice  ((3ow,  1984). 

The  dendrites  themselves  contain  negligible  amounts  of  urea.  As  the  ice 
grows,  we  see  the  formation  of  long  columnar  crystals,  with  each  crystal 
composed  of  several  plate dke  dendrites.  The  resulting  crystal  substructure  of 
ice  plates  and  urea  lamellae  is  in  every  way  analogous  to  the  ice  plate  and 
brine  pocket  substructure  of  the  saline  ice. 


Groove 


Figure  2.1:  Sdiemetic  ol  interface  geometry  of  two  adjacent  crystals  of  sea  ice, 
showing  vertical  dendrite  and  groove  structure  and  horizontal  c* 
axes.  Brine  entrapment  is  confined  to  the  grooves.  The  larger  arrow 
on  the  left  indicates  the  downward  direction  of  growth. 
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The  two  ice  types  app>ear  to  be  very  similar  morphologically  (See  Figure 
2.2).  For  both  ice  types  we  notice  multi-plate  composition  of  individual 
crystals,  liquid  inclusions  along  the  boumdaries  between  the  plates  and  the 
cellular  structure  of  the  individual  plates.  (Gow,  1984) 

Seeded  urea-ice  sheets  characteristically  possess  a  multi-layered  structure 
that  is  best  observed  in  both  vertical  and  horizontal  thin  sections  (See  Figure 
2.2).  A  thin  section  is  a  sheet  of  ice  less  than  1  mm  thick,  mounted  on  a  glass 
slide.  The  top  layer  of  the  ice  sheet  is  called  the  seed  layer  (designated  S),  and 
consists  of  seed  crystals  generally  less  than  1mm  in  diameter  and  up  to  3  cei  ^ 
thick.  Underneath  the  seed  layer  are  one  or  more  layers  of  columnar  crystals 
in  which  the  c-axe$  are  either  vertical  or  horizontal.  Layers  of  mixed 
orientation  are  not  uncommon,  although  individual  layers  tend  to  be 
dominated  by  crystals  of  one  orientation  or  the  other.  The  coltunnar  crystals 
in  this  1  ■  er  known  as  the  incubation  layer  and  designated  ic,  exhibit  crystal 
growth  controlled  by  planar  interface  geometry.  Textures  and  orientations 
encountered  in  the  ic  layer  are  those  typically  observed  in  freshwater  ice. 
Individual  crystals  in  the  ic  layer  contain  very  few  liquid  or  gasesous 
inclusions;  the  resultant  structure  is  optically  transparent,  in  contrast  to  the 
milky  semi-opaque  appearance  of  the  underlying  ice  (Gow,  1984).  By  qmcKiy 
reducing  the  temperature  to  10**F  at  the  time  of  seeding  the  urea  ice  sheet,  we 
can  confine  the  incubation  layer  to  thidmesses  ranging  fix>m  05  to  1  cm. 
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The  incubation  layer  is  und^lain  by  the  dendritic  colunmar  layer, 
designated  dc,  which  is  composed  of  large,  columnar  crystals  of  dendritic 
origin.  This  layer  forms  the  bulk  of  the  thickness  of  the  ice  sheet.  It  consists 
mainly  of  crystals  with  horizontal  c-axes  since  any  crystals  with  c-axis 
orientations  other  than  horizontal  are  usually  eliminated  by  the  process 
known  as  geometric  selection.  This  process  is  well  documented  for  sea  ice 
(Perey  and  Potmder,  1958;  Weeks  and  Lee,  1958). 

The  tips  of  the  dendrites,  which  protrude  down  into  the  solution,  form 
the  final  layer  known  as  the  dendritic  interface.  As  growth  occurs,  substantial 
am.  kjits  of  urea  are  systematically  incorporated  as  layers  of  inclusions 
between  the  ice  plates  in  individual  crystals.  It  is  assumed  that  this  interface 
plays  a  significant  role  in  the  transmission  and  absorption  of  acoustic  energy 
into  the  ice  sheet  The  ragged,  feathery  appearance  of  crjrstals  from  the  bottom 
of  the  ice  sheet  is  evident  in  the  photomicrograph.  (See  Figure  23)  Rigorous 
quantitative  measurements  to  determine  the  exact  length  of  fiiese  tips  have 
not  been  made,  but  a  qualitative  estimation  shows  that  in  saline  ice  these 
doidritic  tips  are  elongated  and  pointed  while  the  dendritic  tips  of  urea  ice 
appear  shorter  and  sturdier  (See  Figures  23, 2.4). 


Figure  2.3:  Photograph  of  basal  thick  section  of  urea  ice  sheet  grown  in 

October,  1989  showing  structure  of  dendrites  and  interdendritic 
grooves.  Amplification  is  X  7. 
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Figure  2.4:  Photograph  of  basal  thick  section  of  saline  ice  sheet  grown  in  1985 
showing  structure  of  dendrites  and  interd^dritic  grooves. 
Amplification  is  X  7.  The  thin  white  lines  are  the  interdendritic 
grooves  while  the  darker  sections  show  the  dendrites  of  pure  ice. 
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Z2  Acoustic  Theoiy: 

2^1  Tnuismission  of  Power- 

Assume  that  the  source  is  a  small  sphere.  It  expands  and  contracts 
uniformly  and  the  signal  travels  outward  radially.  At  normal  incidence,  in 
the  far  field  (klUl)  of  the  circular  piston  transducers  of  the  type  we  were 
using,  the  signal  might  be  assumed  to  be  only  a  function  of  the  distance  R. 

Following  Clay  and  Medwin  (p79-82),  the  wave  equation  for  radially 
propagating  waves  in  a  spherical  co-ordinates  S3rstem  is 


where  p<R,t)  is  die  travelling  pressure  signal.  Guessing  the  answer,  we 
write 

p(R,t)  =  <  b/R)  g(R,t)  (2) 

where  b  is  a  constant 
Substituting  g(R,t)into  (1)  gives 
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8R2 


g(R.t) 


xi_ 

2dt2 


(3) 


which  has  the  same  form  as  the  one  dimensional  plaxw  wave  equation.  Its 
solutions  are  g(R,t) »  g(t-R/c)  and  g(t>R/c). 


We  will  choose  the  outgoing  signal  p(R^t)  » (bc/R)  g(t>  R/c),  in  complex 
form  with  time  dependence  exp(ia)t)  and  write 


p(R,t)  -  (bc/R)  exp(i<o(t-  R/c)) 

Now  we  can  evaluate  b^  by  assuming  that  the  rms  source  pressiue  is  Pq  at 
Rq  ,  therefore  b^  ■  PoRq. 

Alternatively,  the  real  expression  can  be  used 

p(R,t)  ■  (br/R)  cosfojft-  R/c)) 


The  mean  square  value  of  the  cosine  function  is  0.5,  the  source  is  Pq  at  Rq, 
therefore  bf  »  VIPflHo 

Now  we  can  write  the  expression  for  p(R,t)  as 

p(R,t)«  (PqRq/R)  exp(i«(t-  R/c)) 

p(R,t)  -  (VSPqRo/R)  co8(«(t-  R/c))  (4) 


For  a  relation  between  the  aooiistic  pressure  and  the  particle  velocity,  we 
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may  represent  their  relation  in  radial  components 
dp(R,t)  _  ^  du(R.t) 

where  u(R,t)  is  the  radial  component  of  the  particle  velocity.  Substituting 
(4)  into  (5)  and  using  (o  «  kc  yields 


i§!exp{i  (OK-  kR)l  (1+  (6> 

K  IKK 

Since  u  must  have  time  dependence  expficot)  we  use  io>u(R,t) 

dt 

Then  (6)  becomes 

p(R,t)  (1+  «(Pac)  u(R.t)  (7) 

For  large  kR  ( i.e.  many  wavelengths  from  the  source),  u  is  in  phase  with  p 
and  the  relation  is  the  same  for  plane  waves 

p(R,t)  -  (pac)  u(R,t)  kR  »  1 

The  source  transmits  the  power  into  the  water.  We  define  power  as  force 
times  the  velocity  in  the  same  direction.  The  instantaneous  intensity,  the 
power  that  passes  through  a  unit  area,  is  p(R,t)u(R,t)  when  both  quantities  are 
real.  Alternatively  we  use  the  real  part  of  p(R,t)u*(R,t)  when  the  quantities 
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are  complex.  For  the  sound  pressure  Pq  at  Rq  ,  the  complex  pressure  and 
particle  velocity  are 

p(R,t)»  (PqRo/R)  exp(i(o(t-  R/c)) 

u(R,t)  *  (1+  T^^)expU  (£Dt-  kR)l  (8) 

RpAC  ikR  ^ 

The  average  intensity  <I>,  or  the  average  power  passing  through  a  unit 
area^  at  distance  R  is 

<I>  *  Rfi^l  p(R,t)u*(R,t) 

Substitution  of  (8)  followed  by  integration  yields 

8. 

Since  the  beam  pattern  geometry  in  the  near  fidd  of  the  transducers  may 
confuse  the  measurements,  we  take  our  calibration  readings  in  the  far  field 
which  is  defined  as 
R^  2  xa^/X 

by  the  American  Naticmal  Standards  Institute  (Clay  and  Medwin,  p.  155) 
where  a  is  the  radius  of  the  transducer  and  X  is  the  wavelengdt  of  the 


{ 
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transmitting  wave.  At  20  kHz,  the  wavela[\gth  of  the  acoustic  wave  in  water 
is  0.0^  and  the  radius  a  is  0.03m.  Thus,  the  far  field  for  our  circular  piston 
transducer  was  approximately  4  cm.  The  1  /R  calibrations  were  made  with  the 
surface  of  the  water  about  60-80  cm  away  from  the  transducers.  Similarly, 
beam  pattern  measurements  were  made  in  the  far  field  with  the  target 
approximately  50-60  cm  from  the  transducers.  Assuming  that  the  sound 
velocity  in  water  is  1400  m/s,  the  value  of  k  at  20  kHz  is  89.75  m*V  Since  the 
range,  R,  betweoi  the  transmitter  and  receiver  is  about  Im,  we  satisfy  the 
condition  that  the  above  calculations  be  made  where  kR  »  1. 


The  total  power  from  the  source  11  is  the  integral  of  <I>  that  passes 
through  die  surface  of  the  Sphere  of  radius  R 


(10) 


where  dS/R^  is  the  element  of  die  solid  angle. 
Since  I  and  R  are  constant  over  the  spherical  sufiwe 


n«<I>R2|da 


01) 
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This  is  the  power  radiated  through  the  spherical  surface  at  any  radius, 
provided  there  is  no  loss  of  energy  in  the  medium. 

The  rms  pressure  at  the  reference  radius  Rq  is 


In  mks  units,  pressure  is  in  newtcms  per  square  meter  (N/m^) «  pascals 
(Pa),  power  in  watts,  distance  in  meters,  and  px  in  kilograms  per  cubic  meter. 
Generally,  1  m  serves  as  the  reference  Rq  because  it  is  convenient  in 
computation. 

The  rms  pressure  P  at  a  distance  R  is  invmely  proportional  to  the  Ae 
distance  R . 


P» 


(13) 


( Clay  and  hfedwin,  p.  82) 
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2JL2  The  Sonar  equation: 

All  the  characteristics  of  the  sonar  system,  (the  sound  transmission,  the 
scattering  characteristics  of  objects)  are  combined  in  tiie  scmar  equation,  (day 
and  h^win,  p.l03)  After  all  the  quantities  have  been  made  dimensionless 
using  reference  distances  and  referoice  sound  pressure  levels,  the  sonar 
equation  is  usiially  expressed  in  decibels.  Ihus,  the  the  sonar  equation  for 
transmissian  to  a  distant  point  (Le  in  the  far  is 


SPL-SL-TL 

where 

SPL-201og(P/Pr) 

dBrePf 

and 

SL  *201og(Po/Pr) 

dBreP, 

and 

TL  «  20  log(  R/Rq)  +  aR 

where  Pf 

is  the  reference  pressure  (IN/m^),  Rq  is  the  reference  distance 

dm),  a  is  the  attenuation  coefficient  SL  refers  to  the  source  level,  TL  to 
transmission  loss,  and  SPL  to  sound  pressure  level  at  tiie  field  position. 

lfo%vever  when  the  ray  paths  reflect  off  another  surface,  tiie  amplitude  of 
the  reflected  sigiud  is  the  amplitude  of  the  incident  signal  mtiltiplied  by  the 
reflection  coefficient  We  assume  that  the  local  interface  is  planar  with  respect 
to  the  wavdengtii  of  tiie  incident  wave  We  describe  this  reflection  process  in 
the  sonar  equation  by  defiidng  the  bottom  loss  as 
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BL  «  •201ogR|2 


The  net  SPL  is 


SPL-SL-TL-BL 


The  bottom  loss  may  conversdy  be  used  to  compute  the  reflectian 
coefficient  of  the  reflecting  surface. 

We  used  a  system  of  relative  calibration  to  compute  the  reflection 
coefficient  of  the  ice.  In  order  to  refor^ice  the  reflection  coefficients  of  the  ice 
with  a  known  quantity,  we  assume  that  tiie  maeritude  of  the  refl<?ction 
coefficient  off  the  open  water  is  1.  Thus  die  bottom  loss  from  open  water  is  0. 
Now, 


SPLioe  *  SLjq  -TL(q  -BLice  (1) 

SPLwater  *  S^o  water  water  ’®^rater  (2) 

For  similar  distances, 

S^o  ice*  S^o  water 
and  TL|o  TLjq  water 
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Therefore,  subtracting  (1)  from  (2)  we  get 


S^’^'water’  S**Mce*  ®Mce 

20  log  (Pwater/Pr)  -20  log  (Pice/Pr> » ‘20  log  R13 

^ice/^water  *  ^13 

This  method  defines  our  mediod  of  relative  calibration.  By  reading  die 
amplitude  of  the  SPL  off  the  ice,  and  comparing  it  to  die  amplitude  of  the  SPL 
off  open  water  at  a  comparable  distance,  we  can  derive  the  total  reflection 


coefficient,  R13,  for  the  ice. 
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2JL3  The  Plane  Wave  Reverberation  Modeb 

The  reflection  a)effidents  derived  for  the  various  frequencies  as  a 
fimction  of  ice  thidouss  can  be  explained  using  a  well*known  theory  for 
describing  plane  wave  reverberation  by  a  slab  (e.g.  day  atui  Medwin,  p.  64). 

For  a  loss>le8s  slab,  the  ncMrmal  incidence  reflection  ooefricient  is  given  by 

R  P2C2,-P1C1 

P:ft  +  Pici 

where  p  is  die  density  of  the  medium  and  c  is  the  velocity  of  the  acoustic 
wave  in  the  medium,  fri  order  to  introduce  absorption  loss,  we  let  die  sound 
speed  (or  altemativdy  the  wave  number)  in  the  slab  medium  be  compdex  i.e. 
C2«  C2'  iC2”  where  C2’  and  C2’’are  reaL 

Now  consider  an  acoustic  pressure  wave  incident  normally  fn»n  below  on 
a  sheet  of  ice.  The  total  down-travelling  sigzud  is  the  sum  of  the  partial 
transmissions  and  reflections.  Each  padi  within  the  layer  has  a  phase  delay  of 
2Kh  where  K  is  the  wave  number  and  is  defined  as 

K-k-t-ia 
k  m  0f>/c2' 
a»ci)^Vc2’ 
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k  is  the  real  component  of  the  wave  number  and  a  is  the  attenuation 
in  the  medium  in  nepers/m.  The  above  approximations  are  valid  only  when 
I C2'*  I  is  much  less  than  C2'.  By  letting  the  incident  signal  have  unit 
amplitude,  we  can  conclude  that  die  total  reflectimi  is 


R 13  .  Ri2  ♦  Ti2  T2iR23«-^  Z  t  R23R21»‘^ 

Now  assiuning  that 

Ri2  •  -R21 

and  *^12^21*  ^  *  ^12^ 


we  can  express  R13  as 


R.,-  Ri2-KRae-2iKh 
(1+ Ri2  Rjsc-iiKh) 


(Ri2-»-R23e-2i(k4ia)h) 

(l+Ri2R23C-2i(k4ia)h) 
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In  our  case,  the  slab  of  urea  ice  is  bounded  on  the  top  by  a  layer  of  air  and 
on  the  bottom  by  a  layer  of  water.  Assuming  that  there  is  total  reflection  from 
the  top  air-ice  internee,  R23«  *1-  The  magnitude  of  the  reflection  ooeffident 
(R|3)  appropriate  for  the  transducers  in  the  water  is  given  by 

IR,j|«  .  (Itii»  *  e-4ah .  Riie-2«li(2C(»(2kh))) 

(l+R,l*e-<ah .  R,2r2ah(2Cos(2kh))) 

where  R|2  is  the  reflection  coefficient  betweeii  the  ice  and  water,  a  is  the 
attenuation  through  the  ice  in  nepers/m,  h  is  the  thickness  of  the  ice  in  cm 
and  k  is  the  wave-number  in  ice.  Complex  wavenumbers  are  intended  to 
model  absorption  loss.  Because,  the  real  part  of  the  reflection  coefficient  is 
almost  a  factor  of  10  larger  than  the  imaginary  part,  the  magnitude  of  the  real 
part  dominates  die  complex  value.  To  simplify  the  computations,  we  assume 
the  reflection  coefficients  to  be  real.  This  assumption  is  reasonable  for  the 
frequencies  we  are  operating  at  and  the  attenuation  measured  (at  20  Khz,  k«36 
m'^  and  a«3  m'^).  The  model  assumes  uniform  attenuation  throughout  the 
ice  and  no  additional  spreading  loss  owing  to  refraction.  (Jezek  et  aL  1990) 
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2^3.1  LimiUtioiis  of  the  model 

The  plane  wave  reverberation  modd  described  above  assumes  the 
medium  to  consists  of  a  single  layer  widi  uniform  bulk  acoustic  properties. 
However,  previous  experimental  work  (Jezek  et  al.  1988)  deme  on  saline  ice 
has  shown  die  properties  erf  the  ice  to  change  as  a  function  of  depdi.  As  the  ice 
thickens  the  temperature  distribution  within  the  ice  sheet  varies.  This  causes 
changes  in  the  salinity,  porosity  and  permeability  of  die  ice  sheet  Ihe 
structure  and  density  of  the  brine>filled  podeets  show  changes.  Since  Ifae  ice 
rejects  brine  as  die  sheet  gets  thicker,  the  brine  trapped  in  the  early  stages  of 
growdi  of  the  ice  sheet  gets  pushed  out  towards  the  top  of  die  ice  sheet  Since 
brine  is  entrapped  into  the  body  at  the  base,  the  salinity  of  the  sheet  is  highest 
near  the  top  and  bottom  of  a  thick  ice  sheet  Langleben  and  Pounder  (1970)  in 
their  study  of  Arctic  sea  ice  conducted  their  experiments  at  1  m  depth  in  2  m 
thick  ice  and  the  direction  of  propagation  of  the  acoustic  wave  was 
perpendicular  to  the  direction  of  ice  growth.  They  reported  attenuations 
much  lower  than  diose  found  by  Jezek  et  al.  (1969)  for  18  cm  thidc  saline  ice. 
Jezek  et  al.  also  took  normal  incidence  backscatter  measurements  and  found 
that  attenuation  rates  for  6  cm  ice  were  measurably  lower  than  diose  for  3  on 
ice,  substantiating  the  hypothesis  that  attenuatiem  rates  are  diminished  above 
the  dendritic  interface.  Analysis  of  reverberation  data  in  the  8*20kHz  region 
suggests  that  ice  less  dian  a  few  centimeters  thick  is  less  lossy  than  thicker  ice, 
an  observation  that  may  be  associated  with  die  transition  to  a  dendritic 
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interface  (Jezek  et  al  1990).  It  is  recognised  that  the  ice  may  be  modeled  more 
accurately  using  a  multi-layered  model.  The  differoit  layers  may  be  the 
relatively  impervious  incubation  layer,  the  diin  ice  layer,  thicker  ice  layer  and 
the  dendritic  interface. 

Ifowever,  the  ke  sheet  we  worted  widi  was  young  and  relative  thin. 
Because  we  worked  on  young,  thin  ice  we  make  die  assumption  that  the  ice 
sheet  is  a  single,  homogenous  slab.  This  seems  reasonable  given  the 
documented  ice  structure  described  in  section  3.4  where  we  noticethat  the 
dendritic  layer  dominates  the  thickness  of  the  ice.  The  urinity  of  this  layer  is 
estimated  to  be  hdily  constant  at  6  ppt  It  is  also  weakly  supported  by  the  diet 
that  by  using  a  single  layered  plane-vrave  reverberation  model,  we  were  able 
to  obtain  a  reasonably  good  fit,  within  ±  10%,  for  the  theoretical  and 
experimental  data  curves.  Finally,  given  diat  the  ice  sheet  was  grown  only  to 
19  cm  thickness,  the  experimental  data  at  the  lower  frequencies  could  not 
condusitfely  confirm  that  the  attenuation  changed  as  the  ice  grew.  The 
attenuation  parameter  changed  the  rate  of  decay  of  the  reflection  coefficient 
envelope.  The  amplitudes  of  the  peaks  and  troughs  would  determine  the 
slope  of  this  envelope,  and  at  the  lower  frequencies  we  saw  only  2-3  peaks  by 
19  cm  thickness.  A  check  of  the  attenuation  parameters,  showed  that  lowering 
d»  attenuation  parameter  10-15%  provided  a  better  fit  for  the  thinner  ice 
data.  In  addition  the  ice  was  morphologically  relativdy  homogenous  with  no 
obvious  discontinuities,  sudi  as  bubble  layers,  But  would  neccessitate  a 
multi-layered  modeL 
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3.  Experimental  Approach: 

CXir  measurements  took  place  in  December  1989,  at  die  Cold  Redons 
Research  and  Engineering  Labcx’atory's  indoor  cold  pit  in  Hanover,  New 
Hampshire.  The  pit  was  3m  wide  by  3m  long  by  13m  deep.  The  pit  was  filled 
with  tap  water  raised  to  a  urinity  of  10  parts  per  thousand  (ppt)  by  the  addition 
of  finely  granulated  urea.  The  ice  sheet  we  examined  was  the  second  ice  sheet 
grown  for  this  experiment  and  was  seeded  by  spraying  a  fine  layer  of  water 
droplets  evenly  on  the  urea  solution.  The  temperature  in  the  room  %vas 
maintained  between  0^  F  and  10^  for  die  duradon  of  the  experiment  with 
periodic  adjustments  made  to  control  the  rate  of  growth  of  the  ice  sheet  (See 
Figures  3.1, 33).  Consequendy,  the  ice  growth  rate  was  maintained  at  35 
cm/day  for  die  first  10  cm  of  growth  and  then  decreased  to  2-25  cm/day  until 
the  ice  sheet  was  19  cm  thick. 

3.1  Set  up: 

Both  the  growth  and  di^action  experiments  were  performed  using  die 
same  electronic  setup  (See  Figure  35a).  Two  piezodectric  ceramic  transducers 
that  resonated  at  190kHz  with  narrow  beamwidth  (about  9  degrees)  were 
mounted  adjacendy  on  a  horizontal  arm  attadied  to  a  vertical  shaft  The  shaft 
was  then  mounted  cm  a  tripod  aiui  deployed  bdow  die  ice.  The  tripod 
allowed  the  transducers  to  be  rotated  360  degrees. 


Maximum  thickness  19.7  cm 
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Figure  3.1:  Shows  thickness  of  ice  as  a  function  of  time.  Varying  growth  rates 
reflect  changes  in  the  air  temperature  caused  by  power  outages. 


Temperature  Fluctuations  Near  Surface  of  Ice 


T 


were  also  caused  by  power  ouUges. 
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Exp«rimtiital  Sat-up  15,17,20  and  25  kHz 


Function  Ganarator  Powar  Amplifiar  Voltage  transformer 


Oscilloscope  Voltage  amplifier  aP  Filter 


Figure  33:  Etectronics  sel^p  for  die  growth  and  diffraction  phases  of  the 
experiment  at  die  lower  frequencies. 
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Thus,  patches  of  the  iinder-surface  of  the  ice  could  be  izuiependently 
insonified  at  nonnal  incidence.  A  ffewlett  Padcard  (HP)  pulse  gaierator 
8011A  produced  trigger  timing  and  gate  duration  (400  microseconds  for  all 
frequencies)  for  our  pings  and  an  HP  pulse/function  generator  8111A 
produced  tone  bursts  (no.  of  wave<ycles)  at  the  various  resonance  frequencies 
of  our  transducers.  We  used  7  cydes  for  the  lower  frequendes  and  37  cycles  for 
190  kHz.  These  were  die  minimum  number  of  cydes  dut  allowed  die  signals 
to  reach  a  steedy-stats  peak  value.  Hie  output  of  the  function  generator  ( IV 
at  lower  frequendes,  1657  V  at  190  kRr)  was  amplified  through  a  Hsmon' 
Kardon  power  amplifier  that  provided  a  gain  of  30.  Hie  Harmon-Kardon  fed 
a  1:200  voltage  transformer  that  was  connected  to  the  transmitting  transducer. 
Hie  voltage  across  the  transducer  was  measured  to  be  approximately  50n  ^  V. 
The  receiving  transducer  was  connecmd  to  a  AP  Circuit  Corp  Frequency  . uter 
with  a  pass  band  ranging  from  3  kHz  to  24  kHz  for  the  lower  frequendes  and  3 
kHz  to  240  kHz  fcff  die  higher  frequendes.  The  output  from  die  filter  was 
passed  through  an  operational  amplifier  cascade  widi  a  gain  of  2000  (See 
Figure  34 )  befrm  being  read  on  a  Tdctronic  digital  osdllosoope.  The  amplifier 
consisted  of  a  cascade  of  three  741C  operational  amplifiers.  The  first  two  stages 
had  a  gain  of  10  and  die  final  stage  had  a  gain  of  Z  This  provided  a  total 
amplificatum  of  2000  vdiile  still  maintaining  a  bandwidth  of  70  kHz.  A 
Polaroid  camera  was  used  to  {diotograph  the  observed  waveform  on  the 
oscilloscope  screen.  Hie  amplitudes  and  travd  times  of  die  reflected  pulses 
were  manually  recorded. 


Schematic  of  Op-Amp  circuit 
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Figure  3.4:  Schematic  of  the  operational  amplifier  circuit.  Op  amps  had  a 
cumulative  gain  of  2000. 
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Since  the  transducers  had  been  designed  to  resonate  at  190  kHz,  the  setup 
at  that  frequency  was  modified  by  removing  all  the  amplifier  stages  described 
above.  The  transducers  were  directly  connected  to  the  transmit  and  receive 
electrcmics  (See  Hguie  33).  The  received  signal  still  had  a  20-25  dB  tignal  to 
noise  ratio  widiout  ampUficatkm. 


< 
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Experimental  Set-up  190  kHz 


Function  Generator 


Oscilloscope  AP  Filter 


Hgure  33:  Electronics  8efe>up  for  die  growth  phase  of  ttie  e9^>erixnent  at 
190kHz.  Hie  amplifier  drcuitry  has  been  removed  and  die 
transduoos  are  connected  direcdy  to  the  transmit/receive 
electronics. 
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3^  Calibntioiis: 

3^1  lyR  calibntioiit: 

1/R  calibratioiui  were  performed  at  15, 17, 20, 25  and  190  kHz  in  order  to 
verify  that  the  transducers  exhibited  the  predicted  range  dependence  in  the  hu* 
lield.  (Sec  Figmcs  3^3.6e)  At  the  lower  frequencies,  a  7  cyde  normai- 
inddence  signal  (37  cycles  at  190  kHz)  was  reflected  otf  open  water.  Ihe  peak 
amplitude  and  travel  times  for  the  reflected  signals  were  measured  for 
various  distances,  and  the  results  plotted  rni  a  log-log  graph  to  verify  the  1/R 
dependence  of  amplitude  vs  distance.  (See  Table  1)  The  narrow  beuawidth 
(about  9  degrees)  of  the  transducers  operated  at  190  kHz  affected  the  1/R 
calibrations.  Since  the  centres  of  the  transmit  and  receive  transducers  had  a 
space  of  about  7  cm  between  them,  at  dose  range  a  large  part  of  their  beam 
patterns  would  not  overlap.  This  angular  effect  may  have  led  to  a  diminshed 
dependence  of  received  power  on  range  (See  Figure  3.6e).  We  find  a  standard 
deviation  of  23%  from  a  1/R  dependence  at  190  kHz. 

By  carefully  measuring  the  depth  of  die  transducers  below  the  water  and 
then  dividing  by  the  observed  travd  time  of  the  received  signal,  the  speed  of 
sound  in  water  was  computed  to  be  1400  m/s.  The  speed  of  sound  in  the  ice 
was  measured  to  be  approxirrutdy  3800  m/s  (15%).  A  short  pulse  of  3  bursts 
at  190  kHz  was  transmitted  from  imder  die  ice.  The  difterence  in  travd  times 


Hgure  3.6a:  Graph  showing  1/R  range  dependence  of  the  signal  amplitude 
reflected  from  open  water  at  15  kHz.  Data  points  include 
measurements  taken  before  and  after  experiment.  Percentage  error 
-7% 


15  kHz  Pre/Pott  experimentl/R  calft)ration 
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■  Pre  Experknent  □  Pott  Experiment 
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Figure  3.6b:  Graph  showing  1/R  range  dependence  of  the  signal  amplitude 
reflected  from  open  water  at  17  kHz.  Data  points  include 
measurements  taken  before  and  after  experiment.  Percentage  error 
»7% 


17  kHz  Pre/Post  experimentl/R  calibration 


tt  (ms) 


■  Pre  Experiment  □  Post  Experiment 
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.  Igure  3.6c:  Graph  showmg  1/R  range  dependence  of  the  signal  amplitude 
reflected  from  open  water  at  20  kHz.  Data  points  indude 
measturements  taken  before  and  after  experiment.  Percentage  enor 


=  10% 


20  kHz  Pre/Post  experimantl/R  calibration 


tt  (ms) 


■  Pre  Experiment  ^  Post  Experiment 


Figure  3.6d:  Graph  showing  1/R  range  dependence  of  the  signal  amplitude 
reflected  horn  open  water  at  25  kHz.  Data  points  include 
measurements  taken  before  and  after  eqperiment  Percentage  error 
»11% 


25  kHz  Pra/Poat  axparimontl/R  calibration 
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Figure  3.6e:  Graph  showing  1/R  range  dependence  of  the  signal  amplitude 

reflected  from  open  water  at  190  kHz  before  experiment.  Percentage 
error  »  23% 


190  kHz  pre  experiment  1/R  calttvation 


39 


for  the  received  pulses  from  the  water-ice  and  ice-air  interfaces  was  measured. 
The  speed  wais  computed  by  dividing  the  thickness  by  this  di^rence.  The 
change  in  acoustic  velocity  leads  to  refraction  in  die  ice  resulting  in  beam- 
spread  losses  of  approximatdy  2  dB  in  20  cm  thick  ice  at  20  kHz  (Jezek,  verbal 
communi'  .tion).  This  loss,  however,  has  been  ignored  in  the  plane-wave 
reverberation  model. 

After  the  experiment  had  been  conduded,  1/R  calibrations  were  repeated 
to  verify  the  stability  of  the  system  diroughout  the  course  of  the  experiment 
(See  Figufc  3.6a<3.6d).  UnfcHtunately,  the  transducers  had  to  be  removed  from 
the  water  prior  to  the  recaUbratian,  aiui  the  subsequent  re-soaping  and  re¬ 
immersion  increased  the  signal  level  by  -1%  15  kHz,  4%  at  17  kHz;  1S%  at  20 
kHz  and  22%  at  25  kHz. 


Assuming  a  1/R  dependence  of  the  amplitude,  we  calculated  the  product 
of  amplitude  and  travd  times  for  the  different  ffequendes.  A  percentage  error 
was  computed  by  dividing  the  standard  deviation  by  die  mean  value  (See 
Table  3.1).  The  percentage  errors  were  as  follows: 


15kHz 

IZkfis 

20kHz 

25kHz 

ISSMk 

Pre  calibration  4% 

4% 

4% 

4% 

23% 

Post  calibration  10% 

9% 

7% 

5% 

Pre-Pbst  7% 

7% 

10% 

11% 

23% 
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3^  J  Beam  pattern  calibxatioii: 

Beam  Pattern  measurements  vrere  made  to  determine  the  angular 
dependence  of  die  reflected  signal's  amplitude.  The  transducers  vume 
sulnnerged  in  die  pool  and  a  vacuumed  bulb  attatched  to  a  dun  rod  ivas 
suspended  approodmatdy  50>55  cm  direedy  above  die  transducers.  Hie 
transmitting  transducer  was  excited  at  20  kttr.  &i  order  to  minimiee  the 
pulsewidth  of  die  reflected  signal  and  sdll  provide  an  acceptable  signal  levd, 
two  cycles  were  transmitted.  This  provided  a  pulsewiddi  of  about  75>100 
microseconds  while  generating  a  S/N  ratio  of  about  15*20  dB.  Two  sets  of 
measurements  were  takeiL  For  die  first  set  the  bulb  vns  moved  paralld  to  the 
shaft  in  increments  of  5  cm  until  the  reflected  signal  was  no  higher  than  the 
noise  level  of  about  40  mV  (See  ngures  3.7).  For  the  second  set  the  bulb  was 
moved  perpendicular  to  the  shaft  along  the  axis  connecting  the  two 
transducers;  again  in  increments  of  5  cm  until  the  signal  was  indiscernible 
from  the  noise  (See  Fignics  3.8).  The  target  was  placed  in  the  Add  of  the 
transducers  because  die  beam  pattern  geometry  is  complicated  close  to  die 
transducers.  The  results  were  pltMted  on  a  graph  with  gain  against  lateral 
distance  By  knowing  die  verdcal  and  latmal  distances  from  die  transducers  to 
the  target,  we  could  calculate  the  angular  dependence  of  the  reflected  power. 


Beam  pattern  Perpendicular  to  shaft  (pre/post  experiment) 
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Figure  3.8:  Beam  pattern  calibration  showing  target  strength  as  a  function  of 
lateral  distance  parallel  to  die  shaft  away  horn  the  transducers. 
Calibration  was  performed  after  experiment.  The  3dB  loss  shows  the 
half  power  beamwidth. 
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33JL1  Beam  pattem  theory: 


F6r  a  directional  transducer. 


P.£q^  expa(o»-kR) 


where  D  ■  D(8),  die  gaussian  beam  pattern. 

Intensity  away  from  die  source,  at  distance  Rf  from  die  transmitter  is 


Pac  Rt^pac 


Now,  die  power  intercepted  by  die  target  of  cross-sectional  area,  o,  is 


RTPac 


If  we  assume  ^hcrkal  spreadii^  of  the  power  scattered  from  the  target, 
then  the  scattered  power  intensity  is 


n«  — o 

4icR^Pac 


and  die  power  intercepted  by  die  receiver,  at  distance  R^  hrom  die  target 


where  a«3  an. 


Now,  we  have  two  transducers  that  have  their  central  axes  7cm  apart. 

Consider  the  beamwidths  to  be  Gausaan  and  of  die  form  where  Oq 

is  the  half'power  beamwidth  derived  aheee  for  WoAsiacar. 

Now  when  0«6q/2  Le  haif*power  an^e»  then 


.HMeo/20o)2.„ 

.•^'♦.05 
d-  177 


Assume  that  uie  transducers  are  placed  7  cm  apart  and  the  target  is  moved 
52  cm  above  the  axis  connecting  dieir  centres.  The  angle  0^  is  the  angle  made 
between  the  target  and  die  normal  to  the  transmitting  tranducer,  while  angle 

Of  is  the  angle  made  between  the  target  and  the  normal  to  the  receiving 
tranducer.  Now  the  Gaussian  beam  pattern  for  both  transducers  is  of  the  form 

.-d(e,/6o)^.-<J(9,/Oo)2 

*  w 


which  gives  us  the  theoretical  beam  pattern. 
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The  effective  half<power  beamwidth  of  d\e  nuun  lobe  was  the  angle  at 
which  reflected  power  was  half  the  peak  power  i.e.  the  3dB  loss  level  This 
was  calculated  to  be  about  49”  at  20  kHz  (See  Table  3.2).  The  results  of  the 
b«un  pattern  tests  are  summarised  below 


Pre  experiment: 

Measured  8 

Theoretical  0 

Parallel  to  shaft 

50” 

50® 

Perpendicular  to  shaft 

57” 

52® 

Post  Experiment: 

Parallel  to  shaft 

45” 

50” 

Perpendicular  to  shaft 

45” 

52” 

The  half-power  beamwidth  derived  hrom  the  theoretical  plots  is 
approximately  SO'SZ”.  (See  Appendix  2) 

The  lower  measured  beamwiddt  may  be  attributed  to  the  fact  that  the  bulb 
was  not  a  perfecdy  spherical  target  and  its  reflectivity  may  vary  as  the 
function  of  its  positkm  above  die  transducers.  Another  source  of  error  that 
might  account  for  the  skew  in  the  beam  pattern  data  might  be  the  that  the 
transducers  were  not  perfecdy  horizontal  in  their  mounts,  thereby  altering 
the  beampattem. 
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3.4  Obfwrvations  of  Ice  Morphology: 

The  horiz(mtal  and  vertical  thin  section  of  the  ice-sheet  grown  on 
December  13, 1989,  show  no  large  variatioiu  from  the  structure  described 
above.  The  sections  were  prepared  frmn  samples  that  were  removed 
periodically  from  the  same  general  region  on  the  ice  sheet.  Sections  about  0.5 
cm  thick  were  mounted  cm  glass  slides,  and  shaved  using  a  microtome  to  a 
thickness  less  dum  1  mm.  The  crystal  structure  of  the  ice  was  highlighted  by 
photographing  the  secticms  between  crossed  polarizers. 

The  vertical  secticm  shows  features  characteristic  of  urea  ice  as  described  in 
section  2.1.  (See  Figure  3.9).  The  seed  layer  is  barely  distinguishable  at  the  top 
of  the  ice  sheet,  becatise  water  used  to  bond  the  dun  section  onto  the  slide 
melds  the  edges.  The  layor  is  about  0.3-0.4  cm  thick  and  the  layer 
dominates  the  structure  of  thicker  ice.  In  the  case  of  ice  more  dtan  a  few 
centimeters  thick,  we  thereby  assume  the  ice  sheet  to  be  a  uniform  single¬ 
layered  half  space.  We  note  the  Icmg  parallel  lines  of  urea  inclusions  within 
individual  crystals.  There  is  a  progressive  increase  in  die  cross-sectional 
diameter  of  the  columnar  crystals  as  the  ice  thickens.  Ihe  c-axis  orientation  of 
the  crystals  are  generally  widiin  the  horizontal  plane.  The  dendritic  interhioe 
is  not  distinguishable  in  this  section,  since  it  is  very  fragile  and  disintegrates 
during  the  mounting  process. 


A  comparison  of  the  horizontal  cross-sections  shows  how  the  dendrites 
progressively  increase  in  size  with  increasing  ice  thickness.  The  photographs 
taken  in  plain  transmitted  light  shows  vividly  the  urea  inclusions,  air  bubbles 
and  dendrite  boundaries.  The  photographs  taken  between  crossed  polarizers 
reveals  the  different  crystals  by  virtue  of  their  diffmnt  c-axis  orientation  (See 
Figures  3.9-3.12).  The  0.8  cm  thick  ioe  (Sample  A)  shows  the  beginning  of  die 
dendritic  phase  where  the  crystals  are  2-3  plates  diidc.  The  brine  inclusion 
process  has  just  started  and  the  pockets  are  not  wdl  defined.  In  the  3.4  cm 
thick  ice  (Sample  C),  the  dendritic  interface  is  %vell  formed  and  individual 
crystals  are  wdl  defined  and  4-7  dendrites  thidc.  In  the  7.3  cm  thick  ice 
(Sample  E),  the  dendritic  interface  has  sta^lsgd.  and  remains  qualitativdy 
unchanged  for  the  rest  of  the  growth  phase  The  brine  indusions  between  the 
dendritic  plates  are  distinct.  A  fine  line  of  indusions  between  the  plates  is 
indicative  of  a  fine  horizontal  cut  through  the  thickness.  (. 

C  •»  U  c 

A  study  of  tile  plate-widtii  as  a  function  of  thidoiess  revealed  that  periods 
of  high  growth  rate  are  characterised  by  narrow  plate-widtiis  while  periods  of 
slow  growth  are  characterised  by  wider  plate-widtiis  (See  Table  33). 
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Sample 

Thickness(cm) 

Platewidth  (mm) 

Growth  rate  i 

A 

0.8 

0 

33 

B 

14 

035 

33 

C 

3.4 

0.35 

3.5 

D 

5.1 

0.45 

3 

E 

73 

0.55 

1 

F 

8.1 

03 

1 

G 

10.5 

0.45 

3 

H 

116 

03 

0 

I 

14.4 

03 

2 

J 

19.7 

0.55 

1 

Initially,  the  growA  rate  was  3.5  cm/day  and  samples  A,B/  andCshow 
plate-widths  erf  0.35  mm.  As  the  growth  rate  fell  to  2.5cm/day  the  platewidth 
increased  to  0  J  mm.  fin  the  case  of  samples  H  and  J,  where  very  little  growth 
had  taken  place,  the  platewiddi  was  0 J5  mm.  A  diedc  of  the  eff^  of 
platewidth  on  the  reverberation  data  showed  no  observable  oo-relation 
between  the  normal-incidence  reflection  coefficient  and  plate-width. 

Smaller  bubbles  observed  in  the  sections  may  be  air  pockets  entrapped  as 
the  ice  grew  while  the  larger  blobs  (indicated  by  arrows)  are  bubbles  of  air 
trapped  between  die  glass  plate  and  die  ice  during  die  sectioning  process.  In 
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the  19.7an  thidc  ice  (Sample  J)  the  crystals  have  15-20  dendritic  plates.  The 
dendrites  themselves  are  0.55  mm  wide.  No  significant  air  bubble  layors  were 
noticed  and  the  c*axes  orientation  ci  the  crystals  in  die  dendritic  layer  were 
horizontal.  The  brine  pockets,  platelets  and  grains  may  all  act  as  sites  of  stress 
conoentraticm  fiom  passage  of  the  acoustic  pressure  wave  Since  die  crystals 
are  not  perfecdy  vertical  and  do  not  stretdi  the  entire  diickness  of  the  ice 
sheet,  the  aoousdc  wave  may  find  varying  padis  as  it  travels  through  the  ice 
sheet 

The  photomicrograph  of  the  thick  section  of  die  urea  ice  (See  Figure  2J) 
shows  that  the  dendrites  are  0.7mm  thick  and  the  crystals  are  much  larger 
than  the  saline  ice  crystals.  A  quandtadve  measurement  of  the  dimensions  of 
saline  and  urea  ice  dendrites,  grown  under  comparable  conditions,  is  essential 
to  conclusively  prove  diat  the  reflection  coefficients  are  afiected  b  '  ^hese 


structures. 


'igure  3.1BJ:  Samples  A-C  of  urea  ice  sheet.  Sample  A  of  ice  sheet  0.8  cm  thick. 
Sample  B  of  ice  sheet  2  4  cm  thick.  Sample  C  of  ice  sheet  3.4  cm 
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4 .  Growth  Experiment: 

4.1  Procedtuc: 


As  the  ice  sheet  grew,  n<ninal  inddenoe  reflection  coefficient  readings 
were  taken  at  15, 17, 20, 25  and  190  kHz  approximatdy  every  05  cm  of 
thidcness.  The  shaft  was  rotated  to  allow  for  several  independent  areas  to  be 
insonified  for  each  frequency.  The  amplitude  of  the  waveform  was  recorded 
where  the  waveform  showed  a  steady*state,  i.e.  the  amplitude  remained 
tmdumged  regardless  of  the  number  of  cydes  (See  Fig  43a).  The  osdiloscope 
had  two  cursors.  One  cursor  was  kept  constant  throughout  the  experiment  at 
ground  level  to  mark  the  reference  zero  as  well  as  foe  begixuung  of  foe 
transmit  pulse.  The  amplitude  of  the  reflected  signal  was  chosen  by  moving 
the  second  curscv  to  foe  «.eak  of  the  cyde  chosen  to  represent  the  steady  state 
value  (shown  with  arrow).  The  difference  in  foe  peak  value  and  foe  ground 
value  were  read  off  foe  display  on  foe  osdUoscope.  The  travel  time  of  the 
received  waveform  was  recorded  by  moving  the  second  cursw  to  the 
beginning  of  the  first  peak  of  foe  received  waveform.  The  difference  in  travel 
times  between  foe  two  cursors  was  read  off  foe  oscilloscope  display.  Since  foe 
transducers  had  a  beamwidth  of  approximately  40°  at  the  lower  frequendes, 
we  took  4*5  readings  at  each  of  the  lower  fiequendes  and  6>7  readings  at  190 
kHz.  The  mean  value  was  considered  to  be  the  amplitude  of  file  reflected 
signaL  Thidmess  samples  were  taken  at  periodic  intervals  (0.8  cm3>4cm. 
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3.4cm,  5.1cm,  73cm,  8.1cm,  10.5cm,  12.6cm,  14.4cm,  19.7cm)  to  monitar  the 
growth  rate  and  any  changes  in  the  morphology  of  the  ice.  At  the  oondusicr' 
of  the  growth  phase  of  die  experiment,  in  order  to  statistically  validate  die 
accuracy  of  our  reverberation  data,  we  collected  127  distinct  points  of 
reflectkm  coefficient  data  at  190  kHz  and  20  independent  points  of  reflection 
coefficient  data  at  each  of  die  lower  fiequendes.  P.d.f.'s  were  tabulated  for  the 
different  frequendes  and  the  statistical  variations  were  computed  (See  Fig^ues 
4.1a-4.1e). 


Figure  4.1a:  (Cstograms  compiled  the  normalized  echo  amplitude  data 
collected  at  die  end  of  the  gro%^  phase  of  the  expriment  Data  were 
collected  at  15  kHz  when  die  ice  was  19  cm  thick. 

Mean  value  «  0.62 
St.  dev.  a  0.014 
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Figure  4.1b:  Histograms  compiled  from  the  normalized  echo  amplitude  data 
collected  at  die  end  oi  the  growth  phase  the  exprimoit  Data  were 
collected  at  17  kHz  whm  the  ice  was  19  cm  thick. 

Mean  value  «  0.62 
St.  dev. «  0.021 
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Figure  4.1c  Histograms  compiled  from  the  normalized  echo  amplitude  data 
collected  at  die  end  of  the  growth  phase  of  the  expriment  Data  wore 
collected  at  20  kHz  when  the  ice  was  19  cm  thick. 

Mean  value  «  0.16 
St.  dev. »  0.014 
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Figure  4. Id:  Histograms  compiled  from  the  normalized  echo  amplitude  data 
collected  at  the  aid  of  the  growth  phase  of  the  expriment  Data  woe 
collected  at  25  kHz  when  the  ice  was  19  cm  thick. 

Mean  value  >  0.56 
St.  dev.  >  0.012 
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Figure  4.1e:  Histograms  compiled  from  the  normalized  echo  amplitude  data 
collected  at  the  end  of  the  growth  phase  of  the  expriment  Data  were 
collected  at  190  kHz  when  the  ice  was  19  cm  thick. 

Mean  value  «  0.12 
St  dev.  *  0.004 
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4^  Precautions: 

The  surface  of  the  transducer  was  soaped  prior  to  its  immo^sion  in  the 
water  to  prevent  the  formatum  of  bubbles.  The  bubbles  would  create  an 
impedence  mismatch  at  d\e  transducer-water  interface  that  could  advmely 
aH^  the  coupling  of  the  acoustic  signal  into  the  water.  In  order  to  prevent 
the  downward  transmission  of  acoustic  energy  as  well  as  diminish  die  large 
acoustic  coupling  between  the  transducers,  they  were  mounted  on  a  baffling 
that  consisted  of  a  1cm  layer  of  cork,  atop  a  5cm  layer  of  styrofoam.  The  cork 
had  two  circular  cavities  which  snugly  housed  the  transducers.  A  sponge  was 
stuffed  around  die  transducers  to  further  reduce  any  coupling.  To  alleviate 
electronic  coupling,  all  exposed  wiring  was  physically  separated  and  the 
receiving  and  transmitting  electronics  were  placed  apart.  The  coaxial  cables 
from  the  transducers  were  taped  to  the  shaft  to  prevent  entanglement  in  the 
growing  ice  sheet.  As  a  further  precaution  against  cutting  the  wires  while 
freeing  the  shaft  bxnn  the  ice  sheet,  a  45cm  long  circular  sleeve  was  fit  snugly 
around  die  shaft  and  extended  from  15cm  above  the  water  level  to  30cm 
below.  (Sec  Figure  4^ 

When  determining  the  amplitude  of  the  received  pulse,  care  was  taken  to 
sdect  die  right  position  on  the  reflected  waveform.  It  is  impOTtant  to  choose 
that  part  of  die  waveform  diat  represents  the  total  steady-state  reflected  signal 
(See  ^;are  43a).  At  the  lower  frequencies,  since  attenuation  through  the  ice 
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was  low  (0.3  dB/cm),  there  was  considerable  interference  from  the  reflection 
off  the  top  layer.  Fbr  diin  ice,  the  two  reflections  interfered  quickly  and  a 
steady-state  representative  waveform  was  observed  very  early.  Ifowever,  as 
the  ice  thickened,  null  points  appeared  in  the  total  reflection  where  the 
reflections  off  the  top  and  bottom  layer  were  180  degrees  out  of  phase  4See 
Figure  43b).  In  some  cases,  7  cycles  were  not  enough  to  distincdy  see  the 
interfmnce  and  so  the  cydes  were  increased  to  10-12  and  the  amplitude  was 
recorded  where  the  signal  levded  out.  Sometimes,  the  reflected  signal  was 
highly  dependent  on  hrequency  and  small  variations  in  frequency  of  1  mr  2 
kHz  would  alter  the  reflected  signal  by  10-15  dB.  This  occured  infrequently  at 
the  lower  frequencies  but  often  at  ti»  higher  frrequendes. 
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Figure  4.2:  The  physical  setup  of  the  growth  phase  of  the  experiment.  The 
transducers  were  mounted  on  an  arm  susp^ded  from  the  tripod. 
The  electronics  were  housed  outside  the  cold  room  and  connected 
to  the  transducers  using  20  ft  long  coaxial  cables. 


Growth  Experiment  Set  up 
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Figure  4.3a:  Shows  typical  20  kHz  waveform.  Travel  time  measured  from 
beginning  of  pulse,  while  amplitude  was  measured  by  placing 
cursor  at  steady  state  peak  of  waveform. 
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Figure  4.3b:  Sometimes,  due  to  interference  of  the  reflections,  the  steady  state 
value  could  not  be  determined  within  7  cycles.  Hence  the  number 
of  cycles  was  increased  till  the  signal  levelled  out.  Photo  1:  7  cycles 
Photo  2: 10  cycles 
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Figure  4.4  Shows  the  null  caused  by  the  interference  of  reflections  fron\  the 
top  and  bottom  surfaces  of  the  ice  sheet.  Amplitude  of  received 
signal  is  taken  at  the  trough  in  the  signal  after  confirming  its  steady 
state. 
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43  Observations: 

The  results  of  our  growth  experiment  were  rather  surprising.  Considering 
the  similar  morphologies  of  urea  and  saline  ice  and  assuming  that  the 
acoustic  properties  are  independent  of  the  composition  of  the  entrafyiad 
brine,  we  expected  comparable  bulk  properties.  However,  we  observed 
significant  differences  at  the  lower  frequencies. 

As  the  ice  sheet  grew,  normal-incidence  sonar  echo  amplitudes  -were 
recorded  approximately  every  0.5  cm  of  growth.  Using  the  technique 
described  in  the  sonar  equation  theory,  the  total  normal-incidence  reflection 
coefildent  was  derived.  Figures  4.4a-4.4e  show  the  amplitude  of  normal 
incidence  reflection  coefficients  as  a  function  of  ice  thickness  at  difierent 
frequencies.  The  experimental  data  were  fitted  using  a  plane-wave 
reverberation  model  for  a  single-layored  lossy  half-space  boimd  at  the  top  by 
air  and  the  bottom  by  water  (See  Appendix  3).  We  assiuned  total  refiection 
from  the  top  air-ice  interface.  The  model  had  four  variable  parameteis: 
frequency,  attenuation,  bottom-surface  reflection  coefficient  and  velocity  of 
the  acoustic  wave  in  tiie  ice.  The  frequency  parameter  influenced  the  wave 
number  and  thereby  affected  the  number  of  oscillations  observed.  The 
reflection  coefficient  fiom  the  bottom  layer  determined  the  final  steady  state 
value  of  the  total  reflection  coefficient  after  the  oscillations  had  settied  down. 
The  attenuation  parameter  determined  the  exponential  decay  of  the  refiection 
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coefficient  envelope.  The  velocity  parameter  determined  the  spacing  between 
the  troughs  and  peaks. 

As  shown  in  Figure  4.4a,  we  observed  that  at  190  kHz,  the  total  reflection 
for  thin  ice  was  dominated  by  the  reflection  fiom  the  air>ice  interface.  As  the 
ice  thickened,  attenuation  caused  loss  of  signal  fiom  the  the  top  interface.  At 
approximately  8  cm  thickness,  almost  all  the  reflected  signal  was  frmn  the 
bottom  water-ice  interface.  The  limited  data  for  the  thin  ice  allowed  several 
velocities  to  fit  the  data  equally  well.  Since  the  reflection  coefficient  dropped 
off  quite  rapidly  we  assumed  a  high  attenuation  of  2.86  dB/cm.  Since  after 
about  8  cm  thickness,  the  reflection  coefficient  acheived  a  steady  state,  the 
bottom  layer  reflection  coefficient  was  determined  to  be  0.12.  By  measuring 
the  difference  in  echo  arrival-times  from  the  top  and  bottom  layers  for  a  3 
cycle  signal,  the  velocity  was  found  to  be  3800  m/s. 

The  estimated  standard  deviation,  o',  of  the  reflection  coefficient  is 
computed  assiuning  that  it  has  a  normal  distribution. 


where  o  is  the  standard  deviaticm  of  a  particular  data  set  t  is  the 
adjustment  factor  derived  fiom  the  t  distribution  fiiat  gives  us  a  value  that 
the  data  point  has  a  90%  probability  of  lying  within  o'  of  the  mean.  The  value 
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of 't*  is  a  function  of  the  number  of  observations  and  the  probability  that  data 
point  lies  within  the  specified  range.  Values  for  t  were  detennined  using  the 
table  on  p692  of  Probability  and  Statistics  by  Morris  H.  DeGroot. 

On  average  we  took  about  5  measurements  per  thidcness  at  190  itife  The 
variance  of  each  daU  set  was  computed.  The  average  standard  deviation  for 
the  data  sets  was  the  square  root  oi  the  average  variance.  The  adjusted 
standard  deviation  of  a  typical  data  point,  using  die  t  distribution  voe  then 
0.04.  Using  information  from  the  probability  density  function  (p.d.f.),  for  the 
reflection  coefficient  we  computed  a  standard  deviation  of  0.004. 

For  the  20kHz  data  die  attenuation  was  almost  a  factor  of  10  lower  at 
0.31dB/cm  (See  Figure  4.4b).  This  low  attenuation  allowed  for  a  significant 
contribution  of  reflected  signal  from  the  air>ioe  interfiice.  Phase  cancellation 
between  the  upper  and  lower  reflections  caused  a  sharp  drop  in  the  total 
reflection  coefficient  at  about  10cm.  Unlike  widi  the  190  kHz  data  the 
reflection  coefficient  did  not  drop  off  to  a  steady-state  value.  Rather,  we  noted 
an  unusually  high  bottom-layer  reflection  coefficient  of  0.35,  very  comparable 
to  that  of  hard  freshwater  ice  which  has  a  reflection  coefficient  of  0.41.  We 
took  about  3-4  data  points  per  thickness.  The  average  standard  deviation  for 
the  data  points  adjusted  using  the  t-distribution  is  0.08.  The  standard 
deviation  of  die  reflection  coefficient  from  the  p.di  data  was  0.014.  It  might  be 
noted  that  since  the  p.d.f.  data  for  the  lower  frequencies  is  limited  to  20 


70 


samples,  it  may  not  show  the  shape  of  a  nonnal  distribution.  The  velocity 
was  estimated  to  be  3800  m/s  which  gave  us  the  best  fit  to  the  data. 

Results  at  the  other  lower  frequencies  showed  omsistently  hi^  reflection 
coefficients  (See  Figures  4.4c4.4e).  At  15kHz  we  observed  a  reflection 
coefficient  of  0.6,  an  attenuation  of  0.27  dB/cm  and  velocity  of  3800  m/s.  From 
the  p.d.f.  data  we  found  the  standard  deviation  of  the  reflection  coefficient 
was  0.014.  The  average  s.d.  of  a  data  point  was  0.09. 

At  17  kHz  we  observed  a  reflection  coefficient  of  0.4,  an  attenuation  of  0.28 
dB/cm  and  velocity  of  3800  m/s.  The  standard  deviation  of  the  reflection 
coefficient  from  the  p.d.f.  was  0.0214.  The  average  s.d.  of  a  data  point  was  0.08. 

At  25  kHz  wb'  observed  a  reflection  coefficient  of  0.4,  an  attenuation  of  035 
dB/cm.  The  velocity  that  provided  a  best  fit  for  the  data  was  3950  m/s.  The 
standard  deviation  of  the  reflection  coefficient  from  the  p.d.f.  was  0.012.  The 
average  s.d.  of  a  data  point  was  0.1. 
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Figure  4.4a:  Mean  value  of  normalized  ed»o  amplitudes  at  190  kHz  as  a 
function  of  ice  thickness.  Smooth  curve  was  calculated  using  a 
simple  plane  wave  reverberation  model  using  the  parameters 


Attenuation  >  0.31  dB/cm  R.C.s  0.35  Velocity^  3800  m/s 


r 


72 


Figure  4.4b:  Mean  value  of  normalized  echo  amplitudes  at  20  kHz  as  a 
function  of  ice  thickness.  Smooth  curve  was  calculated  using  a 
simple  plane  wave  reverberation  model  using  the  parameters 
shown.  Also  shows  range  of  attenuation  values. 


Attenuation  s  0.28  dB/cm  R.C.s  0.4  Velocity^  3800  m/s 
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Figure  4.4d:  Mean  value  of  normalized  echo  amplitudes  at  17  kHz  as  a 
function  of  ice  thickness.  Smooth  curve  was  calculated  using  a 
simple  plane  wave  reverberation  model  using  the  parameters 
shown.  Also  shows  range  of  attenuation  values. 
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TfliMc  4tl 

Summary  of  bulk  acoustic  values  of  the  ice  at  various  frequencies: 


Freq. 

LSx 

Avg.  S.d. 

Vdocitv 

15  kHz 

0.6 

0.09 

0.27  dB/cm 

3800  m/s 

17  kHz 

0.4 

0.08 

0.28  dB/cm 

3800  m/s 

20  kHz 

035 

0.08 

0.31  dB/cm 

3800 m/s 

25  kHz 

0.4 

0.1 

0.35  dB/cm 

3960 m/s 

190  kHz 

0.12 

0.04 

2.86  dB/cm 

3800m/s 
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4.4  Analysis 

A  comparison  of  experimentally  derived  values  for  urea  ice  showed 
marked  differences  from  diose  derived  from  saline  ice  (See  Figure  4.5a). 
France^  et.  aL  (1988)  found  low  reflection  coefficients  for  cold  saline  ice.  This 
might  indicate  good  coupling  with  the  dendritic  layer.  As  the  ice  warmed  up, 
however,  it  experienced  an  upward  dumge  in  reflection  coefficient  It  is 
speculated  dut  warming  temperatures  may  erode  the  complicated  inters 
connections  between  brine  pockets  formed  during  fieezing,  allowing  freer 
flow  of  brine  between  the  dendrites.  This  might  reduce  the  attenuation 
sufiered  by  the  wave  traveling  duough  the  ice.  Hence,  energy  might  be 
reflected  from  the  stratigraphic  discontinuities  like  bubble  layerSr  bigh&  in 
the  ice  and  ultimately  contribute  to  increased  received  signal  levels.  ( Jezek  et 
al.l989) 

For  example,  Francois  found  duit  at  20  kHz  cold  ice  had  a  reflection 
coefficient  of  0.07.  Warm  ice,  howe>^,  had  a  reflecticm  coefficient  of  0.23,  a 
tripling  in  value.  At  lower  frequencies  our  data  on  urea  ice  showed  even 
higher  reflection  coefficients,  almost  comparable  to  that  of  freshwater  ice's 
refiection  coefficient  of  0.4. 


In  order  to  determine  direedy  die  r^lection  coefficient  from  die  derived 
data,  we  analysed  the  waveform  of  die  received  pulse  at  a  thickness  of  19  cm. 
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The  total  distance  travelled  by  die  sound  wave  in  the  ice  was  thus  38  cm  and 
at  3800  m/s  it  would  take  100  ^s  b^ore  the  reflection  hrom  the  top  surface 
would  interfere  with  the  reflection  ftom  the  bottom  surface.  The  pulse 
received  hrom  the  19  cm  thick  ice  and  die  calibration  pulse  were  digit^ed.  The 
amplitude  of  the  ice-reflected  waveform  was  thra  compensated  for  losses 
based  on  its  arrival  time.  By  comparing  the  peak  values  before  100  |is  of  diis 
waveform  to  the  peak  values  before  100  ps  in  the  calibration  waveforms  we 
derived  the  reflection  coeff  ients  for  die  bottom  layer.  These  values  were 
determined  to  be 
Iit2l£A2 


Frequency 

Digitised  r,c: 

Moddr.c. 

Error 

15  kHz 

0.55 

.6 

-8.3% 

17  kHz 

0.45 

.4 

+12.5% 

20  kHz 

0.4 

.35 

+14% 

25  kHz 

0.45 

.4 

+12.5% 

190  kHz 

0.12 

.12 

As  shown  above,  these  values  agreed  widi  the  values  that  were  used  to  tit 
the  modd  widi  the  experimental  data  within  an  error  of  14  %,  well  within 
the  allowable  range  of  reflection  a)efticients. 


We  do  not  have  a  theory  on  the  exact  mechanism  of  reflection,  but  we 
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propose  that  dues  might  lie  in  the  dendritic  layer  of  the  ice-water  interface 
and  the  coupling  effect  it  has  on  the  transmission  of  acoustic  energy  into  the 
ice.  The  porous  nature  of  this  layer,  higher  urinity,  and  permeability  offer 
good  acoustic  coupling  into  the  ice.  Although  not  quantitatively  measured,  it 
appears  that  the  dendritic  tips  of  saline  ice  and  urea  ice  are  morphologically 
different  The  saline  ice  dendrites  have  tips  that  are  elongated  and  thin  while 
the  dendritic  tips  of  urea  ice  appear  shorter  and  sturdier.  (Sec  Figures  23, 2.4) 
Since  the  dendritic  interface  plays  a  major  role  in  the  incorporation  of  brine 
up  into  the  body  of  the  ice  sheet,  this  might  also  have  implications  on  the 
mechanism  of  brine  incorporation  into  the  two  ice  types.  This  might  lead  to 
changes  in  the  volume  of  urea  entrapped  into  the  ice  sheet.  Coi\sequently,  the 
presence  of  a  layer  of  high  urea  concentration  dose  to  the  dendritic  interface 
may  lead  to  higher  reflected  energy.  Unfortunately,  unlike  salinity 
measurements,  no  convement  method  exists  to  accurately  determine  the 
urinity  of  an  ice  sheet  as  a  function  of  its  thickness. 

Bogorodsky  et  al.  (1975)  proposed  that  most  of  the  attenuation  in  young 
sea  ice  results  from  the  viscous  losses  due  to  the  presence  of  salt  water  in  the 
inter-dendritic  grooves.  The  results  of  experiments  on  elastic  wave 
attenuaticm  in  ice  structures  show  that  the  frequency  dependence  on 
attenuation  is  not  very  high,  although  absolute  attenuation  levels  are  very 
high.  It  is  deduced  that  attenuation  in  sea  ice  is  caused  by  losses  due  to 
viscosity  or  thermal  conductivity,  but  not  by  scattering.  (Bogorodsky  et  al. 
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1975)  Thus  they  determined  that  calculation  of  the  attenuation  coeffidait  in 
sea  ice  gives  the  dependence 

a(dB/m)-7(f(kHz))0-5 

A  range  of  attenuation  values  were  made  to  fit  the  data  (See  Figures  4.4a> 
4.4e).  The  central  value  that  best  Bt  the  overall  data  yielded  a  close  fit  with 
Bogorodsky's  et  al.  predicted  values  at  the  lower  fiequendes.  However,  at 
higher  frequendes  Bogorodsky's  attenuation  was  too  low.  It  has  been 
suggested  by  Jezek  et  al.  (1990)  that  the  changing  structure  of  brine  pockets  due 
to  the  warming  of  the  ice  may  affect  die  attenuation.  A  change  in  the  brine 
incorporation  mechanism  as  a  result  of  morphological  difierences  in  the 
dendritic  interface  may  lead  to  a  difference  in  the  viscosity  o:  ihe  liquid-filled 
pockets.  This  might  explain  the  discrepancy  we  observed  in  attenuation 
levels. 

At  lower  fiequendes,  our  data  were  below  Jezek  et  al.  saline  ice  data. 

Saline  ice  reverberation  data  in  the  8-20  kHz  region  suggests  ice  less  than  a 
few  centimeters  thick  is  less  lossy  than  thicker  ice,  an  observation  that  may  be 
assodated  with  a  transition  to  a  dendritic  interfiioe  that  occurs  after  2  cm 
growth.  ( Jezdc  et  al ,  1989)  This  low  attenuation  allows  considerable  sigiud  to 
be  transmitted  through  the  ice,  contributing  to  very  high  total  reflection 
coeffident  of  almost  0.9  for  the  thin  ice  thicknesses. 
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Over  all  frequency  regimes,  we  measured  attenuation  greater  than 
Langleben  and  Pounder  (1970).  Their  measurements,  however,  wese  made  at 
Im  depth  on  ice  that  was  2  m  thick  and  the  direction  of  propogation  was 
perpjendicular  to  the  growth  direction.  This  is  also  the  low  salinity  .fortion  of 
the  ice  column.  Given  the  high  urinity,  porous  and  permeable  nature  of  the 
dendritic  interface,  we  hypothesize  that  the  dendritic  interface  may  play  an 
important  role  in  the  acoustic  attenuation  process. 
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In  ordor  to  determine  the  value  of  the  attenuation  independent  of  the 
models,  the  calibration  and  19  cm  thick  ice  digitised  waveforms  were  used. 
The  calibration  waveform  was  reduced  by  multiplying  it  with  the  bottom 
surface  reflection  coefficient.  This  gave  the  waveform  of  the  first  reflection  off 
the  bottom  water>ioe  interface.  TKen,  this  waveform  was  subtracted  fix>m  the 
19  cm  thick  ice  waveform.  The  resultant  waveform  was  reflection  off  the  top 
surface  of  the  ice  and  the  subsequent  multiple  reflections. 

Using  the  relationships  of  the  pressure  reflection  coefficients  derived  in 
the  theory  of  the  plane  wave  reverberation  model,  we  have 


Ri2  -  -R21 

and  "^12*^21*  ^  “  ^12^ 

Assume  a  lossly  slab,  bound  at  the  top  by  medium  3  and  at  the  bottom  by 
medium  1.  A  normally  incident  wave  of  pressure  is  incident  on  the 
bottom  surface.  The  1st  reflection  from  the  bottom  surface  is  simply 
Pl-PiRl2 


The  second  reflection  i.e  off  the  top  surface  is 
P2  *  Pj  T|2R23^21®^"*^^ 


where  h  is  die  thickness  of  the  slab  and  a  is  the  attenuaticm. 
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The  third  reflection  is 

P3»  PiTj2R21^23^^21®*P^’*^^ 


Using  these  relationships  we  calculated  the  amplitudes  of  the  inflections 
and  noticed  that  diey  added  up  to  the  total  amplitude  of  die  observed 
waveform.  The  attenuation  values  were  those  used  in  the  model. 


Table  4.3 

Frequency 

Ei2 

2nd  refl. 

15  kHz 

.55 

125mV 

17kHz 

0.45 

380 

20kHz 

0.4 

508 

25kHz 

0.45 

160 

3rd  refl. 

Total 

Observed  2nd» 

22mV 

147mV 

125mV 

55 

435 

420 

72 

580 

700 

12 

172 

220 

The  values  of  the  attenuation  used  in  the  model  are  within  ±  10%  of  the 
values  found  in  the  observed  waveform  data.  We  find  that  the  attenuation 
values  used  in  the  model  for  15  and  17  kHz  are  lower  than  those  derived 
^m  the  waveform.  The  opposite  is  true  for  20  and  25  kHz. 

The  model  values  of  attenuation  seem  to  fit  the  data  of  the  thicker  ice 
better.  By  lowering  the  attenuation,  about  10-20%,  the  model  seems  to  have  a 
better  fit  on  the  diin  ice  section  of  the  data.  This  substantiates  the  hypodiesis 
that  attenuation  rates  change  as  the  ice  grows  thideer. 


15  KHz  2od  *  roWocUons 


17kHz  2n(i*  reflections 


Fig  4.6b:  Shows  the  resultant  waveform  of  the  2nd  and  3rd  reflections  from 
the  ice  sheet  at  17  kHz.  Gives  an  independent  measure  of 
attenuation,  velocity  and  r.c.  R|2  =0.45,Tj2=0.55,T2i=1.45Incident 
signal  =  1.56  V 


20  kHz  2nd  *  reflections 
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Fig  4.6c:  Shows  the  resultant  waveform  of  the  2nd  and  3rd  reflections  from 
the  ice  sheet  at  20  kHz.  Gives  an  independent  measure  of 
attenuation,  velocity  and  r.c.  R12  =0.4,Tj2=0.6,T2i=1.4  Incident 
signal  =  2.36V 
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5.  Diffnction  experiment: 

5.1  DiffracdonTlieoiy 

Diffraction  can  be  qualitatively  imderstood  by  graphical  constructions 
using  Huygen's  principle.  However  in  cxder  to  get  a  quantitative  measuire  of 
the  amplitudes  of  diffracted  waves  we  can  quantify  Huygens  principle  by 
using  the  Hrixnholtz-Kirchoff  (HK)  integral.  The  H-K  integral  relates  the 
wave  field  U  on  the  scattering  surface  to  the  tidd  U(Q)  at  a  point  Q.  We  find 
that  U  at  the  surface  is  approximately  RUj  where  is  the  incident  wave  tidd 
and  R  is  the  reflection  oo^dent  The  H>K  integral  is 


U(Q) 


dn 


dS 


(1) 


where  e'^/r  is  the  "point  source”  or  Green's  function  for  the  Huygen's 
wavdets,  r  is  the  distance  from  dS  to  Q,  d/dn  is  the  derivative  along  the 
normal  to  the  surface  and  the  d/dn( )  is  evaluated  at  dS.  This  equation 
assumes  that  the  source  and  recdver  are  in  the  same  medium  and  that  the 
reflection  coefficient  R  is  constant  over  the  entire  rough  surface  Since  U5 
has  the  form  of  an  expanding  wavetiont  it  has  dimensions  of  m'^and  is 
generally  a  function  of  frequency.  (Clay  and  Medwin,  p  320) 


( 
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Let, 


U,- 


er^ 

r 


(2) 


Then, 

To  find  the  field  strength  at  Q,  we  convolve  U(Q)  with  P(0  ,  the  emission 
spectrum  of  the  source,  at  range  rQ  and  get 


p(t,Q)  ■  ro 


P(0U(Q)ei2*ftdf 


(3) 


For  an  impulsive  source  the  frequency  spectrum  P(0  is  a  constant,  Pg,  wd 
we  get 


p(t,Q)-n)Ps 


I 


U(Q)e^df 


Defining, 


u(a»)-  I  U(Q)e“»*»df 


(4) 
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and  substituting  we  find  that 

p(Q/t) »  ro  Ps  u(Q,t)  (5) 

Since  we  will  only  be  calculating  relative  amplitudes,  we  need  only 
calculate  u(Q).  By  substituting  the  H-K  integral  into  (4)  and  integrating  we 
Hnd  that 


where  the  first  term  is  the  reflected  wave  when  the  source  is  under  the 
surface.  The  second  term  represents  the  difiracted  wave  where  the  H-K 
integral  has  been  transformed  to  a  line  integral  over  the  polar  azimuth,  6, 
measured  in  the  plane  of  the  diffracting  surface.  This  integral  appears  often 
and  we  can  define  D(t)  as  follows  to  obtain  Trorey's  (1970)  form: 


We  can  calculate  the  diffracted  boundary  wave  for  a  simple  plate,  and  then 
construct  more  complicated  surfiues  by  joining  a  number  of  plates  together 
and  superimposing  the  solutions  of  the  individual  plates  (Jezek,  1980).  Let  the 
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source  be  on  the  vertical  z  axis  and  a  flat  plate  of  width  r|'  r2  meastired  on  the 
x-axis,  where  the  range  vectors  R|  and  R2  are  measured  from  the  source  at  Q 
to  the  edges  of  the  plate.  The  plate  extends  to  plus  and  minus  infinity  on  the 
y-axis.  The  difi^acted  wave  from  the  edge  at  X|  is  geiterated  as  R|  moves  from 

9s  0  to  ic/2  and  simultaneously  as  R|  moves  from  6*  0  to  -ic/2.  Thus,  fix-  the 
edge  we  have 


T^  is  the  time  corresponding  to  the  horizontal  range  to  the  difiractor 
divided  by  wave  velocity  and  T|  is  the  range  to  the  diffractor  divided  by  the 
wave  velocity.  As  t  tends  to  T| ,  D(t)  becomes  infinite  and  has  an  impulsive 
response.  In  order  to  estimate  the  amplitude  of  arrivals  when  tsT|  we 
average  the  integral  over  a  small  time  interval  A  .  Letting  t*  T|  -t-A  and 
simplifying  by  dropping  the  second  order  terms,  we  find  that 
for  T|  St  STi+ A 
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CKt)  -  2(Ti2  A)r'  tan*i  (8) 

T* 

Equation  (7)  is  the  solution  for  a  semi-infinite  strip  when  the  source  is 
over  the  strip  and  the  other  edge  is  very  far  away.  In  our  experiment,  this  is 
the  case  when  the  transducers  are  under  the  ice  sheet. 

For  the  case  where  the  source  is  not  over  the  strip,  and  both  edges  are 
illuminated  then 

u(Q)«B^(Di(t)-D2(t)) 

IK* 

e  die  specular  reflection  is  missing,  and  the  integral  over  8  has  two 
boundary  waves.  The  boundary  wave  from  the  near  edge  of  the  strip  (R|)  has 
a  positive  sign,  and  the  boundary  wave  from  the  far  edge  (R2)  has  a  negative 
sign. 
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5^  Procedure 

For  the  diffraction  phase  a  slot  was  cut  in  the  ice  when  it  was  19  cm  thick. 
The  slot  was  1.84  m  long  and  13  m  wide  Particular  eftbrt  was  spent  on 
making  sure  that  the  edges  of  the  ice  were  cut  as  strai^t  as  possible.  The 
transducer  shaft  was  removed  from  dte  tripod  and  mounted  on  a  sliding 
apparatus  that  allowed  die  transducers  to  be  moved  gradually  from 
completely  under  the  ice,  to  under  the  edge,  to  completely  under  the  open 
water.  The  transducers  were  exdted  at  20  kHz  and  short  bursts  of  200  ^s  were 
transmitted.  The  pulse  generator  put  out  a  signal  of  2  cydes. 

To  ensure  that  the  transducers  were  horizontal,  the  transducer  arm  was 
mounted  such  that  it  followed  a  3.6°  slope  as  it  traversed  the  length  of  the 
sliding  apparatus.  Accordingly,  compensation  had  to  be  made  to  adjust  for  the 
difference  in  depdi  as  the  transducers  were  moved.  The  total  change  in  depth 
for  the  entire  lengdi  of  the  transducer  traverse  was  8.5cm.  (Sec  F^ure  5.1a,b) 

Measurements  of  the  reflected  and  diffracted  echo  arrivals  were  made 
every  5  cm,  from  85  an  away  from  the  edge  under  die  ice  to  50  cm  away  frmn 
the  edge  under  the  open  water.  The  first  waveform  was  termed  >85  cm,  the 
waveform  directly  under  the  edge  was  termed  Oom,  while  the  last  waveform 
was  termed  50  cm  (Sec  Figure  5.1c). 
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Figure  5.1a:  Di^action  experiment  set  up  showing  the  position  of  the 

transducers  and  their  motion  under  the  edge.  The  transducers  were 
moved  in  increments  of  5  cm. 


DIffraciion  Experiment  Set  Up 


Lateral  distance  versus  Surface  distance 


Lateral  Distance  from  edge  (cm) 
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Figure  5.1c:  The  position  of  the  tranducers  when  completely  under  the  ice 

was  tamed  -85  cm.  The  transducer  was  moved  out  towards  the  edge 
in  increments  of  5  cm.  The  position  when  the  tranducers  were 
under  the  edge  was  termed  0  cm.  When  the  transducars  were 
furthest  to  the  left,  the  position  was  termed  50  cm. 


Diffraction  Experiment  CoH>rdinate  eyetem 

ICE 


180  cm 


Transducer  position 


50  cm 


0  cm 


-85  cm 
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A  nttmber  of  reflections  were  received  and  difft’acted  echo  from  the  edge 
of  the  ice  sheet  was  isolated.  While  the  transducers  were  tmder  die  ice  sheet, 
echoes  were  received  frcnn  the  ice  surface,  the  edge  of  the  ice  sheet,  die  edge  of 
the  water,  and  probably  from  the  sidewalls.  While  the  transducers  were  under 
the  water,  echoes  were  received  from  the  edge  of  the  ice,  the  edge  of  the  water, 
the  lateral  edges  and  probably  the  sidewalls.  (Sec  Appendix  4).  The  travel 
times  of  reflected  die  diffracted  arrivals  were  recorded  by  moving  die  second 
cursor  to  the  beginning  of  die  wawef(»m  of  interest 

Since  the  pulsewidth  was  almost  2(X)  microseconds,  the  direction  pulse 
received  from  the  edge  merged  with  the  reflection  from  the  ice  sheet  when 
the  transducers  were  about  25  cm  away  from  the  edge,  under  the  ioe.  When 
the  transducers  were  under  the  watn  the  diffracted  arrival  could  be  detected 
until  about  20  cm  from  the  edge.  Beyond  that  range,  its  arrival  could  be  noted 
but  it  was  swamped  by  the  reflection  from  the  water  before  the  signal  peaked 
(SceFigfue5.2). 

The  difrracted  arrival  from  the  water  edge  was  also  noticed  from  about  -45 
cm  until  the  edge.  Once  the  transducers  were  underwater,  the  diffracted 
arrival  from  the  water  edge  was  swamped  by  the  reflection  from  open  water 
(Sec  ^nre  S3a,b).  The  reflections  from  the  surface  of  tiie  ice  and  the  surface 
of  tile  water  were  easily  observed,  because  of  tiieir  rdativdy  large 
amplitudes.  (See  Rgiirc  5w4a,b) 
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In  order  to  isolate  the  directed  ptilse  from  the  total  waveform  we 
digitised  the  waveforms.  The  diffracted  waveforms  were  digitised  using  the 
TEK  4010  digitiser.  Software  written  by  Lisa  Belec  of  Thayer  School  of 
Engineering  was  used  to  download  the  data.  The  data  was  manipulated  and 
plotted  using  Microsoft  Excel  2.2 

The  waveforms  were  assumed  to  be  linear  between  peaks  and  so  only  the 
peaks  were  recorded.  The  -35  cm  reflection  data,  (which  was  taken  35  cm  away 
from  the  edge  imder  the  ice)  was  considered  to  be  purely  reflection  from  the 
ice.  This  waveform  was  then  superimposed  and  subtracted  from  the  other 
(-35  to  0  cm)  waveforms.  The  resultant  waveform  was  assumed  to  be  the 
diffracted  arrival  from  the  edge  of  tlw  ice.  Similarly,  the  50  cm  waveform  (50 
cm  away  from  the  edge,  under  the  water)  was  considered  to  be  excliisively 
reflection  from  water.  This  waveform  was  subtracted  from  die  5cm  to  45  cm 
waveforms,  to  isolate  the  diffracted  arrival. 

There  were  several  constraining  factors  that  limited  the  sophistication  of 
this  approach.  First,  die  photographs  of  die  waveform  was  enlarged  optically 
using  a  photocopier.  The  resulting  enlargement  was  then  mounted  on  a  grid 
that  determined  the  coordinates  of  the  pmnts  duit  were  to  be  digitised. 

Hnally,  the  waveform  was  manually  traced  by  diddng  die  cursor  on  the 
appropriate  peaks  and  troughs.  Each  stage  of  diis  transfer  of  data  might  result 
in  errors  in  the  amplitude  and  travel  times  of  the  signal  waveform. 
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Figiure  5.2:  The  diffracted  signal  was  observed  directly,  without  interference 
only  when  the  transducers  w^e  vnder  water.  Under  the  ice,  the 
reflection  from  the  ice  interfered  with  the  directed  arrival 


Diffraction  Travei  times  for  ice  edge 
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Figure  53a:  The  diffracted  signal  hx>m  the  water  edge  was  observed  directly, 
without  interference  only  when  the  transducers  were  under  ice. 
Under  the  water,  the  reflection  from  the  open  water  interfered  with 
the  diffracted  arrival 


Diffraction  Travei  times  for  water  edge 
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Figure  5.3b:  The  amplitude  of  the  signal  from  the  water  edge  could  be 

observed  directly  when  the  transducers  were  under  the  ice.  The 
amplitude  increased  gradually  as  we  got  nearer  to  the  edge, 
qu^tatively  following  the  expected  pattern. 


Diffraction  from  water  edge 
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The  next  stage  of  the  data  reduction  process  involved  isolating  the 
diffracted  arrival  by  subtracting  the  waveforms  from  each  other.  The  *35  cm 
waveform  was  considered  to  consist  of  only  an  ice  reflection.  This  waveform 
was  subracted  from  the  -30  to  0  cm  waveforms  by  aligning  the  Hrst  peak  of  the 
two  signals  (Sec  Figuze  5  J).  The  following  peaks  were  then  coi\sidered 
aligned  because  no  abrupt  chaz^  in  phase  was  observed,  and  any 
discrepancies  were  attributed  to  the  data  transf^  process.  The  signals  were 
then  subtracted  and  the  resultant  waveform  was  considered  to  be  the 
diffracted  arrival  since  there  were  no  other  arrivals  expected  in  that  time 
frame. 

For  the  case  where  the  transducm  were  below  the  water  the  50  cm 
waveform  was  considered  to  comprise  of  only  an  open-water  reflection.  This 
waveform  was  subtracted  horn  the  5-30  cm  waveforms.  The  diffracted  arrival 
was  observed  before  the  reflection  foom  die  water.  Hence,  the  two  waveforms 
could  no  longer  be  aligned  by  comparing  their  Hirst  peaks.  Instead,  thdr 
maximum  peaks  were  aligned,  assuming  that  the  reflection  Horn  the  water 
was  steady  enough  and  large  enough  relative  to  the  diffracted  arrival  to  be 
unaffected  by  the  interference  of  the  diffracted  arrival  (See  Figuze  5.6).  The 
most  significant  limitation  of  this  approach  was  that  by  pre-aligning  the 
peaks,  phase  dianges  in  the  di^cted  signal  arrival  could  not  be  accounted 
for.  In  addition,  because  only  die  peak  values  had  been  recorded,  the 
subtracticm  process  was  extremely  sensitive  to  the  proper  alignmoit  of  the 
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peaks.  In  effect,  a  45**  phase  difference  could  dramatically  alter  the  results  of 
the  subtraction  of  two  waveforms;  h’om  a  subtraction  between  two  peaks  to  a 
subtraction  of  a  trough  from  a  peak. 


Figure5.5:  Shows  subtracted  waveform. 

a) Only  ice  surface  reflectkm  waveform  (*35  cm) 

b)  under  ice  waveform 

c)  resultant  difiracted  arrival 


(a)  *35  cm 


Figure5.6:  Shows  subtracted  waveform. 

a) Only  opoi  water  reflection  waveform  (50  cm) 

b)  under  water  waveform 

c)  resultant  diffracted  arrival 
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5 J  Obsenratioiia: 

The  amplitude  of  the  di^cted  signal  was  recorded  by  visual  observation 
of  the  the  received  waveform.  However,  the  pulsewidth  of  the  signsl  (almost 
200  ps),  caused  the  reflection  from  the  the  ice  sheet  to  dominate  d\e  diffracted 
arrival  when  the  transducers  got  within  25  cm  laterally  from  die  edge.  In 
addition,  interference  bet%veen  the  arrivals  from  die  water  edge  andihe  ice 
edge  was  observed.  When  the  transducers  were  under  die  %vater  the 
beginning  of  the  diffracted  arrival  was  observed  upto  about  20  cm  away  from 
the  edge,  after  whidi  most  of  the  difracted  arrival  was  swamped  by  die 
reflection  off  the  open>water.  Qose  to  die  edge  there  was  considoahie 
interference  between  the  reflection  from  the  open-water  and  die  edgs  of  the 
water  as  well  as  the  bottom  of  the  ice  sheet 

In  order  to  iscdate  die  diffracted  arrival  from  the  wavefcmn,  we  employed 

% 

the  subtraction  procedure  described  above.  When  die  transducers  were  under 
the  ice,  the  diffracted  signal  was  recovered  by  subtracting  the  under-ice  signsl 
from  the  total  reflected  signal.  Similarly  when  the  transducers  were 
completely  under  water,  die  diffracted  signal  was  recovered  frcHn  the  total 
signal  by  careful  subtraction  of  the  t^wn-water  signal  from  the  total  reflected 
signal 

In  mder  to  cmnpare  these  values  we  devdoped  a  diffraction  modd  on  the 
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basis  of  the  Uieory  in  Section  5.1.  For  the  case  where  the  transducers  were 
under  the  ice,  it  was  assumed  that  the  ice  sheet  was  an  infinite  strip  with  one 
edge  being  the  diffiacting  edge  and  the  other  edge  being  at  infinite  distance. 
Assuming  the  incident  signal  to  be  an  impulse  function,  the  values  for  the 
reflection  from  the  ice  sheet  were  derived  using  the  first  part  of  equation  7. 

To  calculate  the  peak  value  of  the  impulse  response  tite  directed  signal, 
equation  8  was  used  to  find  the  instantaneous  value  of  D(t»T),  and  then 
combined  in  equation  7  to  derive  the  peak  value  of  the  diffracted  signal.  The 
signal  was  normalized  by  dividing  the  directed  arrival  by  the  reflected 
arrivaL  This  gave  us  a  ratio  that  when  multiplied  by  the  reflectimi  fiom  the 
ice  gave  us  tiie  theoretical  value  of  the  diffracted  signal  Hnally,  the  signal 
was  compensated  for  the  beam  pattern  efiects  of  the  transducers  to  get  the 
theoretical  value  received  at  the  receiver.  These  values  were  compared  to  the 
peak  values  that  we  observed  and  fhe  values  derived  from  the  subtraction. 
(See  Hgoics  5.7a,  5.7b,  5.7c) 

There  are  several  assumptions  made  in  the  model  in  order  to  simplify  the 
ccnnputation.  The  modd  is  strictly  only  valid  for  the  region  where  the 
transducers  are  under  the  ice  sheet.  Once  the  transducers  are  under  the  water 
that  the  diffracting  edge  is  no  longer  illumined  on  only  one  side  In  order  to 
properly  account  for  the  diffraction,  the  vertical  face  of  the  ice  sheet  would 
have  to  be  taken  into  consideration.  This  would  require  that  a  value  for  D(t) 
be  computed  for  this  second  sheet  and  then  superimposed  cm  the  value 
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computed  for  foe  horizontel  sheet.  However,  foe  effective  surface  area  of  the 
vertical  sheet  is  very  small  compared  to  the  total  area  illumined  by  the 
transducers.  We  have  assumed  in  the  model  that  foe  source  is  a  spherically 
radiating  smirce.  A  compensation  for  beam  pattern,  which  assumes  that  the 
entire  edge  is  equally  illumined,  has  been  made  However,  the  beam  pattern 
of  the  transducers,  causes  the  different  positions  along  foe  edge  to  be 
illumined  at  different  intensities.  The  assumption  may  be  justified  because  at 
20  kHz  foe  beam  pattern  is  apprmdmately  50^,  whidi  is  sufficiently  wide  The 
model  assumes  that  foe  source  is  an  impulse  function.  According,  it 
calculates  the  impulse  response  of  foe  reflected  and  diffiacted  arrivals. 
However,  the  amplitude  of  D(t)  at  t»T^  is  infinity.  We  handle  this  infinity 
and  foe  source  function  infinity  as  shown  in  Equation  8.  The  value  of  A,  foe 
time  interval,  affects  foe  computation  of  D(t)  since  0(t)  is  inversely 
prop<ntional  to  foe  square  root  of  A.  The  value  of  foe  reflected  signal  is 
inversely  proportional  to  A.  The  net  effect  is  that  increasing  A  causes  foe 
diffiracted  values  to  increase  at  a  more  gradual  rate.  (Sec  Hguics  5.7a,  5.7b,  5.7c) 
Assuming  A  to  be  foe  pulsewidfo  of  our  signal  (i.e.  200ns)  we  were  able  to  get 
a  reasonably  good  fit  within  ±  10%,  with  foe  subtraction  data.  Surprisingly 
the  observed  data  fit  foe  model  very  well  for  the  region  where  the  transducers 
were  bdow  the  water.  This  might  indicate  that  foe  effect  of  foe  vertical  edge 
does  not  fonn  a  substantial  part  of  foe  total  diffi'acted  arrival. 


Diffraction  Values  delta  «  lO^is 
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Hgure  5.7a:  Thecwetical  values  of  diffracted  arrival  using  a  lime  step  of  lOps 
compared  to  'bo  observed  arrivals  and  the  values  of  the  diffracted 
arrivals  obtained  by  subtraction  analysis. 


Diffraction  Values  daita  -SOtis 


Figure  5.7b:  Theoretical  values  of  diffracted  arrival  using  a  time  step  of  SOps 
cmnpared  to  the  observed  arrivals  and  the  values  of  the  diffracted 
arrivals  obtained  by  subtraction  analysis. 


Diffraction  Values  della  -200^s 


Figure  5.7c:  Theoretical  values  of  diffracted  arrival  using  a  time  step  of  200ps 
ccunpared  to  the  observed  arrivals  and  the  values  of  the  diffracted 
arrivals  obtained  by  subtraction  analysis. 
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6.  Enor  Analysis: 


The  measurement  error  in  the  1/R  calibration  data  could  be  averaged  to  be 
about  10%.  This  error  accounts  for  dw  instability  we  noticed  in  die  system  for 
the  duration  of  the  experiment  In  fact,  considering  diat  the  transducers  had 
been  removed  from  the  water  and  soaped  before  re-calibration,  the  actual 
calibration  error  may  be  even  less.  Ihe  reflection  coefficient  K  is  a  functimi 
of  two  independent  variables 

Signal 

Calibration 

The  signal  errar  and  the  calibration  error  are  considered  statistically 
independent  processes.  Taking  the  derivative,  and  eliininating  the  cross 
ternns,  we  find  that 


%  error* 


dK 

R 


F6rl90kHz, 


negligible.  This  is  determined  by  knowing  the  shape 


of  die  distribution  to  be  Gaussian  and  thereby  determining  the  error  in  the 
meait  The  specific  process  to  derive  the  standard  deviation  is  derived  in  the 
Section  43.  Therefore  the  resultant  error  is  simply. 


d£al. 

Cal 


23% 
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But  this  can  largely  be  attributed  to  the  foct  that  the  narrow  total 
beamwidth  of  the  190  kHz  transducers,  2Sm‘^(1.15/ka)  >  52**,  coupled  with 
the  finite  spacing  betuveen  them  leads  to  a  variation  in  the  1/R  dependence 
when  the  transducers  are  cloee  to  the  water  surface. 

At  the  lower  frequencies  we  note  that  the  reflection  coefficient  does  not 
lend  itself  to  a  percentage  error  calculation.  The  standard  deviations  o£  the 
data  points  are  siuninarised  as  follows: 


Table  6.1 

£isav9hq. 

S.4.(fPf) 

15  kHz 

0.09 

0.014 

17kHz 

0.08 

0.021 

20kHz 

0.08 

0.014 

25kHz 

0.1 

0.012 

190kHz 

0.04 

0.004 

The  errors  in  the  reflection  coefficient  may  arise  from  various  factors  in 
the  system.  The  error  in  measuring  the  o>rrect  travel  time  and  amplitude 
owing  to  misplacing  the  cursor  on  the  waveform,  especially  when  die 
waveform  does  not  show  a  steady  state  sigiud,  might  lead  to  catastrc^hic 
errors.  The  lower  error  in  the  190  kHz  data  may  be  attributed  to  the  fiict  that 
the  transmit/receive  electronics  was  devoid  of  amplifiers  that  have  an  error 
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associated  with  them. 


All  the  reflection  coefficient  measurements  were  pegged  against  the  value 
of  the  reflection  coefficient  off  the  open  water.  Even  though  only  one 
measurement  was  made  at  0  thickness,  the  calibration  value  was  vary  stable. 
The  calibration  errors  for  the  Iowa*  frequencies  are  summarised  below. 

Table  6.2 


Fr«wgvcy 
15  kHz 
17kHz 
20kHz 
25kHz 


Calil 


grorpggentags 

7% 

7% 

10% 

11% 


This  mor  takes  into  account  the  robustness  of  the  system  during  the 
course  of  the  experiment.  It  may  be  biased  high  since  the  transducers  were 
removed  flom  the  water  before  recalibradon,  necessitating  resoaping  of  their 
surflice  This  might  have  led  to  the  hi^er  amplitudes  observed  in  the  post 
experiment  calibration. 

The  plane-wave  reverberation  modd  used  to  predict  the  bulk 
acoustic  properties  ctf  the  ice  sheet  was  sensitive  to  die  various  parameters.  At 
190  kHz,  a  wide  range  of  vdodties,  frmn  3200  m/s  4200  m/s,  fit  die 
theoredcd  plot  to  die  data.  Because  of  the  limited  number  of  data  points  at 
the  thin  ice  thicknesses  (upto  3  cm),  the  attenuation  could  be  varied  between 
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2.5  dB/cm  to  3.4  dB/an.  The  value  of  the  reflection  coefficient  was  relativ^y 
stable  around  1.2  ±  0.025.  One  additional  facet  of  error  that  might  not  be 
reflected  in  the  model  is  the  lateral  shift  of  the  data  owing  to  errors  in  the 
thickness  of  the  ice.  Over  the  entire  pool  the  ice  varied  in  thickness*!  1  cm 
from  the  area  that  we  collected  our  samples  fiom.  As  the  transducers  swept 
an  arc  of  almost  1  m  radius  the  ice  thickness  may  have  varied  from  0.5-1  cm. 
This  might  help  explain  some  of  the  discrepancy  in  the  position  of  the  peaks 
and  troughs  in  the  thin  ice  thickness  data. 

For  the  lower  frequencies,  the  model  was  much  more  sensitive  to  changes 
in  the  velocity  parameter.  The  velocity  could  be  varied  between  3650  m/s  and 
3900  m/s  before  a  noticable  difierence  in  the  model  fit  could  be  observed.  At 
25  kHz  the  vdodty  could  be  varied  between  3900^100  m/s.  The  attenuation 
at  the  lower  frequencies  was  very  low  and  variation  in  the  attenuation 
parameter  caused  die  sharpest  differences  in  the  values  of  the  reflection 
coefficient  at  dte  nulls.  A  range  attenuations  is  used  to  fit  the  modd  to  the 
theory.  At  the  lower  frequencies,  the  attenuation  values  that  fit  die  thin  ice 
data  best  were  about  10'15%  lower  than  those  that  fit  the  thicker  ice  data.  This 
indicates  that  the  ice  could  be  better  modeled  using  a  muld-layered  model. 
However,  the  attenuaticm  change  was  not  very  large  A  summary  of  the 
attenuation  values  used  to  fit  the  data  is  given  below: 
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TaMfc43 


Frequency 

Range  of  attenuation  values 

Model  value 

15  kHz 

0.27-0.43  dB/cm 

0.27  dB/cm 

17kHz 

0.17-0.43  dB/cm 

0.28  dB/cm 

20kHz 

0.21-0.34  dB/cm 

0.31  dB/ cm 

25kHz 

0.26-0.43  dB/cm 

0.35  dB/cm 

The  values  of  the  reflection  coefficients  varied  between  ±0.5  of  their  final 
values.  The  values  of  die  reflection  coefficients  of  the  bottom  layer  derived 
from  the  subtraction  data  varied  within  ±14  %  of  those  derived  using  the 
model. 

For  the  beam  pattern  analysis,  we  observed  a  strong  correlation  between 
the  theoretical  and  measured  beamwidths.  In  the  case  of  the  beam  pattern 
tests  done  parallel  to  the  shaft  on  which  the  transducers  were  attatdted,  we 
noticed  that  the  measured  loss  was  within  1  dB  of  the  theoretical  estimates.  In 
the  case  of  the  beam  pattern  tests  dcme  perpendicular  to  the  shaft,  we  noticed 
that  the  measured  power  values  were  2-3  dB  higher  than  the  theoretical 
predictions  when  the  target  was  above  the  receiving  transducer.  When  the 
target  was  over  die  transmitting  transducer,  the  measured  values  were  2-3  dB 
lower  than  the  theoretical  predictions.  This  might  be  owing  to  the  shape  of 
the  bulb  which  was  not  a  perfecdy  spherical  target  Its  reflectivity  may  have 
varied  as  a  function  of  its  position  over  the  transducers.  Overall,  the 
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difference  in  theoretical  and  measured  beamwidths  was  about  10%. 

Although  die  diffraction  model  was  strictly  valid  only  under  the  ice  sheet, 
the  model  values  fit  the  data  fairly  closely,  especially  when  allowanoes  are 
made  in  the  value  of  the  time  interval,  A.  Changing  the  value  of  A  smoothed 
the  curve  of  die  diffraction  amplitude  as  it  peaked  just  below  the  edge  of  the 
ice.  The  subtraction  method  we  used  was  subject  to  limitations  described 
above,  the  most  severe  of  which  were  experienced  when  the  transdoeen  were 
under  the  water.  Hie  results  of  the  subraction  analysis  are  wifiiin  20<25  %  of 
the  theoretical  predictions  when  the  transducers  are  under  die  ice  sheet. 
However,  under  the  water  they  are  almost  an  order  of  2  lower  than  the 
theoretical  values.  The  observed  valtus  of  the  diffraction  data  follow  the 
opposite  trend.  Since  die  diffraction  signal  was  swamped  by  the  ice  reflection 
25  cm  from  the  edge,  no  direct  observations  were  possible.  Ifowever,  when 
the  transducers  were  below  the  water,  the  diffracted  signal  could  be  observed 
directly.  Using  a  A  of  50  ^s,  the  theoretical  model  predictions  were  Mndtin 
±10%  of  the  observed  value  of  the  diffracted  signal. 
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7.  Conclusion: 

In  conclusion,  we  propose  that  the  dendritic  water-ice  interface  plays  a 
critical  part  in  the  attenuation  and  nomud  incidence  reflection  of  acoustic 
energy  horn  the  ice  sheet.  A  small  variation  in  the  morphology  of  this 
interface  may  cause  a  significant  variation  in  the  bulk  acoustic  properties  of 
the  ice.  A  qualitative  study  of  the  interface  reveals  urea  ice  dendrites  that  are 
much  wider  and  sturdier  than  die  saline  ice  dendrites.  Although  die 
morphology  of  the  urea  and  saline  ice  seem  similar,  urea  ice  has  significandy 
higher  reflecdon  coefficients  than  saline  ice  at  the  tens  of  kHz  range.  In  the 
hundreds  of  kHz  region,  the  data  is  more  comparable  to  contemporary  studies 
done  on  saline  ice.  Attenuation  rates  are  lower  than  those  observed  for  saline 
ice.  We  hypothesize  that  the  difference  may  be  a  result  of  dissimilarities  in 
the  concentration  of  the  liquid  filled  podcets,  due  to  the  differencein  the 
brine  incorporation  mechanism  between  the  two  ice  types. 

From  the  above,  it  might  therefore  be  concluded  that  urea  ice  is  a  viable 
simulant  for  saline  ice  in  acoustic  experiments  at  the  hundreds  of  kHz 
frequencies  but  may  not  be  a  suitaUe  substitute  at  the  tens  of  kHz  frequencies. 
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A  comparison  of  the  normal  incidence  refelection  coefficients  observed  for 
several  types  of  cold  saline  ice  and  freshwater  ice  at  188kHz  is  summarised 
below.  Laboratory  data  were  collected  using  facilities  at  the  US  Army  G>ld 
Regions  Research  and  Engineering  Laboratory.  Held  data  for  saline  ke  types 
were  collected  in  the  Fram  Strait  while  data  on  fresh  lake  ice  were  collected 
on  a  small  lake  near  CRREL.  Total  measurement  error  for  all  data  t^res  was 
±10%. 


TflMe71 

Ice  Type 

Reflection  Coefficient 

Uboratorv 

Fi£ld 

Saline 

New 

0.04 

0.09 

Cellxilar 

0.12 

0.16 

Altered 

0.27 

Slush 

0.42 

Fresh  Lake 

1 

0.4 

Urea 

15kHz 

0.6 

17kHz 

0.4 

20kHz 

035 

25kHz 

0.4 

190kHz 


0.12 
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The  phase  vdodty  of  the  aootistic  wave  in  ice  was  found  to  be  3800m/s  ± 
5%.  This  velocity  is  in  agreement  with  that  observed  by  others  as  summarised 
in  the  table  from  Wallerstedt  et  al.  bdow: 


Table  7.2: 

Buimey 

3840 

Bennet 

3830 

Clarke 

3817 

Bogorodski 

3800 

Roethlisber^rr 

3960 

McCammon 

3593 

&  McDaniel 

3832 

Wallerstedt 

3387 

Vahanvaty 

3800 

Hnally,  the  diffraction  values  were  in  agreement  with  those  predicted 
using  simple  Helmholtz-iOrdtoff  theory.  It  appears  that  die  vertical  edge  of 
an  ice  lead  mn.Y  not  contribute  significandy  to  normally  incident  acoustic 
backscatter  from  die  edge. 
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8.  Further  work: 

Much  further  wcn’k  could  be  done  in  order  to  bettter  interpret  the  results 
of  this  experiment  The  plane-wave  reverberation  model  nuy  be  enhanced  by 
assuming  a  multi-layered  ice-sheet.  This  would  allow  a  better  fit  between  the 
data  and  theoretical  predictions  since  fitere  would  be  a  greater  number  of 
parameters  to  vary.  We  could  then  account  for  the  ice  properties  like 
attenuation  that  change  as  a  functicm  of  thidcness.  In  order  to  make  a  more 
rigorous  comparison  of  the  urea  ice  results  with  saline  ice,  the  experiment 
could  be  repeated  using  the  identical  set-up  and  equipment,  using  a  sheet  of 
saline  ice. 

Tests  to  measure  the  urinity  profile  of  the  ice  sheet  might  determine  infact 
if  there  is  a  difierence  in  the  brine  inclusion  processes  of  urea  and  saline  ice. 

A  more  quantitative  study  of  the  denoritic  interfoce  would  also  substantiate 
the  hypothesis  that  there  is  a  qualitative  difierence  in  the  dimensions  of  the 
dendritic  interface  of  the  two  ice  types. 

Better  transducers,  wifii  restmances  that  aUow  greater  signal  to  noise  ratios 
and  shorter  pulsewidths,  would  greafiy  improve  die  results  of  the  diffiaction 
phase  of  this  experiment  This  would  also  allow  us  to  study  die  change  in 
backscatter  fiom  the  edge  of  a  lead,  as  die  acoustic  properties  of  the  ice  duuige 
as  a  function  of  its  thidcness. 
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Further  work  might  include  developing  a  more  sophisticated  signal 
processing  tedmique  to  isolate  the  dittracted  arrival  from  the  reflected  pulse. 
As  an  improvement  I  would  definitely  recommend  automating  the  data 
gathering  process.  Not  <mly  would  tills  alleviate  the  kmg  and  tedious  task  of 
monitoring  tiie  growing  ice  sheet  but  it  might  also  aUow  digitisaticm  of  tiie 
waveforms,  thereby  permitting  more  sophisticated  signal  processing  of  the 
received  pulses. 


1/RCalibration  table 


Pre  experiment  Calibration 


l5Khz 

ITKhz 

20kHz 

25Khz 

190kH2 

travel  time  p: 

Voltage  (m\ 

(V) 

(V) 

(V) 

(mV) 

840 

800 

2 

3.16 

1.24 

125 

1215 

520 

1.32 

2.2 

0.9 

190 

1145 

560 

1.44 

2.36 

0.94 

200 

1080 

640 

1.5 

2.5 

0.94 

172 

1011 

680 

1.68 

2.64 

1 

200 

945 

690 

1.9 

3 

1.04 

212 

875 

800 

2.04 

3.4 

1.1 

216 

800 

840 

2.2 

3.56 

1.24 

180 

900 

800 

2.16 

3.7 

1.28 

1140 

560 

1.6 

2.6 

1.12 

1078 

620 

1.76 

2.72 

1.08 

1008 

700 

1.52 

3.28 

1.28 

940 

800 

1.68 

3.7 

1.32 

870 

640 

2.08 

3.92 

1.48 

800 

760 

2.24 

4.2 

1.6 

Post  Experiment 

Calibration 

tt  15kHz  17kHz 

20kHz 

25kHz 

190kHz 

Mean  tt 

976.47 

Mean  Pre 

691.25 

1.76 

2.85 

1.05 

Mean  Post 

697.14 

1.86 

3.45 

1.31 

Mean  Prepost 

694.00 

1.81 

3.13 

1.17 

186.88 

Table 
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1  /RCaiibration  table 


Travel  time* 

voltage 

e 

iSKhz 

i7Kh2 

20kHz  25Khz 

190kHz 

Pre  Experiment 

672000 

1680 

2654.4 

1041.6 

105000 

631800 

1603.8 

2673 

1093.5 

230850 

641200 

1648.8 

2702.2 

1076.3 

229000 

691200 

1620 

2700 

1015.2 

165760 

687480 

1698.48 

2669.04 

1011 

202200 

652050 

1795.5 

2835 

982.8 

200340 

700000 

1785 

2975 

962.5 

169000 

672000 

1760 

2848 

992 

144000 

Post  Experiment 

720000 

1944 

3330 

1152 

638400 

1824 

2964 

1276.8 

668360 

1897.28 

2932.16 

1164.24 

705600 

1532.16 

3306.24 

1290.24 

752000 

1579.2 

3478 

1240.8 

556800 

1809.6 

3410.4 

1287.6 

606000 

1792 

3360 

1280 

STOeviations  Pre 

24654.91 

74.27 

115.55 

45.65 

42531.04 

STOeviations  Post 

68206.28 

155.13 

216.82 

59.45 

STOeviations  Pre  Pos 

47985.62 

119.61 

304  59 

124.26 

42531.04 

Mean  Pre 

668466.25 

1698.95 

2757.08 

1021.86 

165768.75 

Mean  Post 

664165.71 

1768.32 

3254.40 

1241.67 

Mean  Prepost 

666459.33 

1731.32 

2989.16 

1124.44 

185768.75 

Error  percentages  - 
Pre 

stdev/  mean 

0.04 

0.04 

0.04 

0.04 

0.23 

Post 

0.10 

0.09 

0.07 

0.05 

Preoost 

0.07 

0.07 

0.10 

0.11 

0.23 

Table 


Beam  pattern  analysis 


Post  calibration  Beam  pattern  analysis 


Perpendicular  to  shaft 

depth- 

52.29 

lat  dist  amplitude: 

tt  ps  tot.  distance 

P/PO 

angle 

power  dB 

Theoretical 

40 

68 

997 

65.8 

0.68 

37.41 

-3.32 

-6.19 

35 

82 

950 

62.9 

0.75 

33.80 

-2.48 

-5.16 

30 

88 

900 

60.3 

0.74 

29.84 

-2.61 

-4.12 

25 

116 

842 

58.0 

0.90 

25.55 

-0.90 

-3.09 

20 

120 

812 

56.0 

0.87 

20.93 

-1.20 

-2.13 

15 

100 

785 

54.4 

0.68 

16.01 

-3.29 

-1.28 

10 

120 

780 

53.2 

0.79 

10.83 

-2,08 

-0.60 

5 

136 

747 

52.5 

0.87 

5.46 

-1.22 

-0.16 

0 

158 

747 

52.3 

1.00 

0.00 

0.00 

0.00 

•5 

124 

777 

52.5 

0.79 

-5.46 

-2.03 

-0.14 

-10 

108 

782 

53.2 

0.71 

-10.83 

-2.99 

-0.57 

-15 

92 

800 

54.4 

0.63 

-16.01 

-4.01 

-1.25 

-20 

88 

827 

56.0 

0.64 

-20.93 

-3.90 

-2.11 

-25 

76 

875 

58.0 

0.59 

-25.55 

-4.57 

-3.08 

-30 

0 

0 

60.3 

0.00 

-29.84 

-4.12 

-35 

0 

0 

62.9 

0.00 

-33.80 

•5.16 

-40 

0 

0 

65.8 

0.00 

-37.41 

-6.19 

Par8''*»l  to 

Shaft 

depth- 

52.29 

f 

52 

912 

62.9 

0.50 

33.80 

■5.51 

•5.69 

iO 

68 

897 

60.3 

0.59 

29.64 

-3.92 

-4.18 

25 

92 

832 

58.0 

0.74 

25.55 

-1.98 

•2.90 

20 

115 

805 

56.0 

0.87 

20.93 

-0.6€ 

-1 .86 

15 

120 

775 

54.4 

0.85 

16.01 

-0.78 

-1  .C4 

10 

144 

762 

53.2 

0.98 

10.83 

0.43 

•0.46 

5 

165 

757 

52.5 

1.10 

5.46 

1.38 

-0.12 

0 

152 

722 

50.5 

50.54  0.93 

0.00 

0.00 

0.00 

-5 

128 

752 

50.8 

0.85 

-5.65 

-0.82 

-0.12 

-10 

112 

756 

51.5 

0.77 

-11.19 

-1.73 

-0.46 

-15 

100 

722 

52.7 

0.72 

-16.53 

-2.31 

-1.04 

-20 

64 

805 

54.4 

0.49 

-21.59 

-5.66 

-1.86 

-25 

52 

825 

56.4 

0.43 

-26.32 

-6.82 

•2.90 

-30 

0 

870 

58.8 

0.00 

-30.69 

-4.18 

-35 

0 

0 

61.5 

0.00 

-34.70 

-5.69 

Baamwidth 
Meaaurad  Theoretical 
Perpendicuiar  45  50 

Parallel  45  52 


■  Jl 

Table 


Beam  pattern  analysis 


Pre  calibration  Beam  pattern  analysis 

Perpendicular  to  shaft 

depth- 

52.5 

lat  distal  amplitude:  tt 

ps  real  distance 

P/PO 

angle 

power  dB 

Theoretical 

40 

72 

985 

66.0 

0.73 

37.30 

-2.74 

-6.19 

35 

80 

930 

63.1 

0.74 

33.69 

•2.61 

-5.16 

30 

92 

877 

60.5 

0.78 

29.74 

-2.13 

-4.12 

25 

128 

820 

58.1 

1.01 

25.46 

0.06 

-3.09 

20 

135 

795 

56.2 

0.99 

20.85 

•0.08 

-2.13 

15 

140 

787 

54.6 

0.97 

15.95 

-0.26 

•1.28 

10 

152 

780 

53.4 

1.01 

10.78 

0.08 

-0.60 

5 

152 

760 

52.7 

0.98 

5.44 

-0.15 

-0.16 

0 

156 

750 

52.5 

1.00 

0.00 

0.00 

0.00 

-5 

152 

777 

52.7 

0.98 

-5.44 

-0.15 

-0.14 

-10 

120 

797 

53.4 

0.80 

-10.78 

-1.97 

-0.57 

-15 

112 

805 

54.6 

0.78 

-15.95 

•2.20 

-1.25 

-20 

96 

835 

56.2 

0.70 

-20.85 

•3.04 

-2.11 

-25 

84 

872 

58.1 

0.66 

•25.46 

-3.60 

•3.08 

-30 

52 

925 

60.5 

0.44 

-29.74 

-7.09 

-4.12 

-35 

56 

950 

63.1 

0.52 

•33.69 

-5.70 

-5.16 

-40 

0 

0 

66.0 

0.00 

•37.30 

-6.19 

Parallel  to 

Shaft 

depth- 

53.2 

35 

0 

0 

63.7 

0.00 

33.34 

•5.69 

30 

76 

895 

61.1 

0.70 

29.42 

-3.15 

-4.18 

25 

88 

825 

58.8 

0.75 

25.17 

-2.54 

•2.90 

20 

100 

812 

56.8 

0.79 

20.60 

•2.02 

-1.86 

15 

124 

798 

55.3 

0.93 

15.75 

•0.63 

-1 .04 

10 

136 

790 

54.1 

0.98 

10.65 

-0.19 

•0.46 

5 

144 

780 

53.4 

1.01 

5.37 

0.08 

-0.12 

0 

144 

760 

53.2 

1.00 

0.00 

0.00 

0.00 

-5 

144 

785 

53.4 

1.01 

-5.37 

0.08 

-0.12 

-10 

128 

802 

54.1 

0.92 

•10.65 

-0.72 

-0.46 

-15 

120 

837 

55.3 

0.90 

-15.75 

-0.92 

-1 .04 

-20 

96 

862 

56.8 

0.76 

-20.60 

-2.37 

-1.86 

-25 

72 

902 

58.8 

0.61 

-25.17 

•4.29 

-2.90 

-30 

60 

937 

61.1 

0.55 

-29.42 

-5.21 

-4.18 

-35 

0 

0 

63.7 

0.00 

•33.34 

-5.69 

Beam  width 
Measured  Theoretical 
Perpendicular  57  50 

Parallel  50  52 


Table 


DIM  h  (500),  r(500) 

OPEN  #1;  name  "moosedat",  access  output,  create  newold,  organization  text 
ERASE  #1 

PRINT  "Alpha  -> 

INPUT  alpha 

PRINT  "Velocity  ->  "; 

INEUT  vice 

PRINT  "Coefficient  ->  "; 

INPUT  rl2 

PRINT  "Frequency  ->  "; 

INPUT  freq 

I£T  k2  =  (2*3.14*freq)/(vicen00) 

FCR  j  =  1  to  220  step  1 
lET  h(j)  =  j/lO 

LET  u  =  rl2^2  +  e3?>(-4*h(j)*alpha) 

LET  u  =  u  -  rl2’*exp(-2*alpha*h(j))*2*cos(2*k2*h(j)) 

LET  d  =  1  +  rl2''2  *  e:?3(-4*h(j) ’^alj^) 

LET  i  =  d  -  rl2*exp(-2*alpha*h(j))*2*cos(2*]c2*h(j)) 

LET  r(j)  =  (u/d)^.5 

SET  #1:  POINTER  EM) 

PRINT  #1:  h(j);chr$(9);r(j) 

NEXT  j 


CLOSE  #1 
EM) 


diffract  theory  values 
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